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1.  INTRODUCTION 


1.1  Background 

The  U.S.  Army  has  an  ongoing  requirement  for  accurate  and  timely  forecasts  of 
low-level  turbulence  (LIT)*  in  support  of  the  operation  of  helicopters  and 
other  low  flying  aircraft,  especially  in  areas  of  complex  terrain.  Some 
recent  reports  of  accidents  and  the  loss  of  flying  time  during  training  exer¬ 
cises  because  of  poor  LIT  predictions  have  defined  a  need  to  evaluate  current 
forecast  procedures  and  to  determine  whether  the  development  of  better  tech¬ 
niques  is  feasible. 

1.2  Objectives 

As  prescribed  in  the  "Statement  of  kiork  (TCN;  87-597),"  the  objectives  of  the 
current  study  are  "...to  assemble  a  data  base  of  concurrently  measured  surface 
and  airborne  turbulent  intensities;  compile  a  listing  of  all  known  forecasting 
methodologies  for  the  prediction  of  mechanical,  thermal,  and  lee  wave  turbu¬ 
lence;  and  utilize  this  existing  information  to  develop  practical  and  user- 
friendly  prognostication  rules  for  LIT  occurrences  that  can  adversely  affect 
low  flying  aircraft. 

1.3  Genera]  Approach 

The  objectives  listed  above  were  accomplished  through  information  from  a 
general  literature  review  and  a  careful  study  of  the  current  LIT  forecast 
problem  at  the  U.S.  Army  National  Training  Center  (NTC)  at  Fort  Irwin, 
California.  The  literature  review  was  performed  to  isolate  current  and  poten¬ 
tial  LLT  forecast  methodologies,  to  locate  available  data  bases  for  the  future 
development  of  improved  statistical  forecast  techniques,  and  to  investigate 
the  applicability  of  artificial  intelligence  (AI)  and  related  systems  to  the 
LLT  forecast  problem.  The  Fort  Irwin  study  involved  an  on-site  problem  eval¬ 
uation,  the  development  of  a  prototype  LLT  fore  cast/ nowcast  system,  and  the 
development  of  a  data  base  to  test  and  further  improve  the  proposed  system. 
Finally,  a  series  of  recommendations  have  been  developed  on  the  basis  of  the 
combined  results  of  the  literature  review  and  the  NTC  study. 


2.  TECHNICAL  DISCUSSION 

2.1  Literature  Review 

2.1.1  Background 

The  dimensions  of  those  atmospheric  motions  that  adversely  affect  aircraft  in 
flight  are  a  function  of  aircraft  design  and  speed.  Critical  response  scales 
commonly  range  from  a  few  tens  to  a  few  hundreds  of  meters.  In  the  atmo¬ 
spheric  boundary  layer,  where  there  is  often  active  turbulent  exchange  of  heat 
and  momentum  between  the  surface  and  the  atmosphere,  the  typical  dimensions  of 
turbulent  eddies  are  proportional  to  the  height  above  the  ground;  that  is. 


*In  this  report  "low-level  turbulence  (LLT)"  is  defined  as  bumpiness  in  flight 
within  the  planetary  boundary  layer. 
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they  lie  in  the  range  of  the  strongest  response  for  most  aircraft.  For  this 
reason,  and  because  of  the  slower  speeds  and  restricted  maneuverability  of  low 
flying  aircraft,  the  prediction  of  LLT  is  one  of  the  most  important  tasks  for 
an  aviation  forecaster  supporting  for  low  level  flight  operations. 

In  the  present  section,  the  literature  has  been  surveyed  to  determine:  (1) 
the  current  state  of  the  art  of  LLT  forecasting,  (2)  the  possibility  of 
developing  improved  LLT  forecast  techniques,  and  (3)  the  applicability  of  AI 
to  the  improvement  of  LLT  techniques. 

2.1.2  Current  LLT  Forecasting  Procedures 

LLT  forecasting  procedures  now  in  use  recognize  the  inability  of  operational 
data  networks  to  resolve  rnesoscale  and  microscale  spatial  and  temporal  charac¬ 
teristics  of  LLf.  Tnerefore,  with  the  exception  of  the  occasional  pilot 
report  (PIREP),  all  current  methods  of  LLT  diagnosis  and  prognosis  deduce  LLT 
from  the  presence  or  some  larger  scale  circulation  which  is  assumed  to  gen¬ 
erate  LLT,  or  from  the  value  of  some  large-scale  parameter,  such  as  a  dimen¬ 
sionless  number  or  index  which  is  related  theoretically  or  empirically  to  LLT 
(for  example,  Burnett,  1970;  Lee  et  al.,  1979;  AWS,  1979;  FAA,  1977;  Mathews, 
1985).  A  ulocx  diagram  of  the  general  LLT  forecast  procedure  is  presented  in 
figure  1. 

The  strong  dependence  of  LLT  forecasts  on  the  prediction  or  observation  of 
specific  larger  scale  circulations  is  emphasized  in  the  literature  by  the 
separation  of  the  iMjority  of  descriptions  of  forecast  procedures  according  to 
the  cause  of  the  LLT  (for  example,  FAA,  1987).  The  major  causes  are  dry 
convection  (thermals),  moist  convection  (thunderstorms,  downbursts,  etc.), 
mechanical  mixing,  mountain  waves,  and  fronts.  Some  of  these  "causes"  occa¬ 
sionally  overlap  or  are  slightly  ambiguous  (for  example,  wind  shear  associated 
with  large-scale  fronts  versus  wind  shear  related  to  thunderstorm  gust 
fronts);  however,  they  are  common  categories  that  appear  throughout  the  lit¬ 
erature  and  will  help  focus  the  discussion  to  follow. 

In  figure  1,  pattern  recognition  generally  refers  to  the  identification  of 
synoptic  weather  patterns  that  support  the  occurrence  of  one  or  more  of  the 
causes  of  LLT  listed  above.  Favorable  large-scale  patterns  for  moist  con¬ 
vective  phenomena  are  wel 1 -described  by  Miller  (1972),  Doswell  (1982,  1985), 
and  Ray  (1986).  Those  patterns  that  support  widespread  dry  convection  are 
discussed  in  detail  in  the  literature  related  to  forecasting  for  gliding  (for 
example,  see  Lindsay  and  Lacy,  1976;  Bradbury  and  Kuettner,  1976;  kiallington, 
196b). 

Synoptic  patterns  that  are  conducive  to  strong  winds  and  mechanical  mixing 
have  been  descrioed  extensively  for  the  continental  United  States,  by  Waters 
(1970).  Patterns  associated  with  mountain  wave  turbulence  are  generally  well- 
known  and  have  been  summarized  by  Alaka  (1960),  Nicholls  (1973),  and  many 
otners.  Large-scale  frontal  patterns  are  also  well-known  from  the  general 
meteorological  literature  (for  example,  Petterssen,  1956;  Palmcn  and  Newton, 
1969;  Keyser,  1986). 
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Figure  1.  Low-level  turbulence  forecast  procedure. 


A  potential  source  for  the  documentation  of  synoptic  patterns  associated  with 
LLT  are  terminal  forecast  manuals  (TFM)  for  base  weather  stations.  Examples 
for  the  Fort  Irwin  area  are  the  NTC  Forecaster  Handbook  (  1987),  Farnham  and 
Gould  (1956),  and  Farnham  and  Vercy  (1969).  A  general  listing  of  TFMs  is 
given  in  AWS  publications  TC-85/001  (1985). 

A  conceptual  model  (figure  1)  is  defined  here  as  a  mental  picture  of  a  meso- 
scale  phenomenon  that  allows  the  forecaster  to  deduce  the  unobserved  LLT  from 
the  wel 1 -observed  larger  scale  pattern.  It  aids  the  forecaster  in  the  inte¬ 
gration  of  sparse  data  into  a  coherent  mesoscale/inicroscale  pattern.  The 
individual  model  is  usually  a  mesoscale  circulation;  it  may  be  based  on  theory 
or  on  an  average  of  special  field  observations,  or  simply  on  experience. 
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Probably  tne  best  example  of  socft  a  model  is  the  thunderstorm  (Palmen  and 
Newton,  19b9;  Atkinson,  1981;  Doswell,  1982,  1985;  Kessler,  1985;  Weisman  and 
Klemp,  198b;  Fujita,  1985)  in  which  LLT  is  associated  with  downbursts,  gust 
fronts,  wind  shears,  and  related  phenomena.  Models  for  dry  convection  are 
discussed  extensively  in  the  soaring  literature  cited  above,  and  by  Scorer 
(1978).  Most  conceptual  mountain  wave  models  used  currently  by  forecasters 
are  a  product  of  tne  Sierra  Wave  Project  (for  example,  see  Alaka,  1960,  for  a 
summary).  More  recent  publications  by  Nicholls  (1973),  Lester  and  Fingerhut 
(1974),  Lilly  (1978),  and  Ourran  (1986a)  have  dealt  with  models  of  mountain 
wave  systems  that  produce  Strong  Downslope  Windstorms  (SOW),  Conceptual 
models  used  by  forecasters  for  deducing  mechanically  induced  LLT  or  for  deter¬ 
mining  the  presence  of  LLT  in  the  vicinity  of  fronts  (aside  from  gust  fronts) 
do  not  nave  a  clear  mesoscale  component. 

In  figure  1,  parameter  evaluation  refers  to  the  process  of  quantifying  the  LLT 
nowcast/forecast  by  determining  the  values  of  critical  parameters.  In  the 
automated  forecast  (for  example,  at  a  weather  central),  this  step  is  accom¬ 
plished  first.  That  is,  once  the  required  data  have  been  acquired  and  ana¬ 
lyzed,  all  "critical  parameters"  nay  be  evaluated  by  computer,  assuming  they 
lend  themselves  to  computation  at  grid  points.  This,  of  course,  produces  a 
nowcast.  A  similar  evaluation  may  be  done  with  predicted  fields.  If  the 
forecast  process  is  manual  (for  example,  at  a  local  forecast  office),  the 
parameter  evaluation  usually  follows  the  pattern  recognition  step  (figure  1). 

Parameters  currently  used  in  LLT  prediction  are  of  three  types:  the  basic 
meteorological  variables,  their  temporal  and  spatial  derivatives,  and  certain 
combinations,  such  as  physical  and/or  empirical  indices,  and  some  measure  of 
terrain  roughness.  Some  of  the  most  common  parameters  are  listed  in  table 
1.  Those  parameters  associated  with  the  prediction  of  LLT  associated  with 
moist  convection  are  extensive  and  well-known,  and  are  not  listed  here.  See 
Miller  (1972)  and  kay  (1906)  for  the  discussion  of  a  wide  variety  of  param¬ 
eters,  indices,  and  other  forecast  tools  (radar  and  satellite  infornetion) 
useful  in  tne  diagnosis  of  LLT  and  wind  shear  associated  with  moist  convec¬ 
tion. 

Forecast  aids  for  LLT  associated  with  large-scale  fronts  generally  depend  on 
some  measure  of  the  intensity  and  speed  of  the  front.  Richwien  (1979)  indi¬ 
cates  that  a  front  with  a  horizontal  temperature  difference  of  at  least  ID  °F 
and  moving  at  30  knots  or  greater  is  associated  with  significant  LLT.  Also, 
it  is  well-known  that  fronts  moving  across  rough  terrain  are  almost  always 
associated  with  LLT. 

Smaller  scale  sea  breeze  fronts  and  convergence  zones  not  associated  with 
thunderstorms  are  known  to  produce  significant  lift  for  gliders  (Wallington, 
1966);  therefore,  they  may  be  a  source  of  LLT  for  some  types  of  aircraft. 
Again,  the  soaring  literature  provides  excellent  guidance  in  the  prediction  of 
those  phenomena  (for  example,  Bradbury  and  Kuettner,  1976). 
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TABLE  1.  COMHOMLY  USED  LLT  FORECAST  PARAMETERS 


DRY  CONVECTION  (THERMAL)  LLT 

Air  temperature 
Temperature  lapse  rate 
Potential  temperature  lapse  rate 
Thermal  Index 
Showalter  Index 
Richardson  Number 


MECHANICAL  LLT 

Surface  wind  speed  and  gusts 
Gradient  level  wind  speed 
Mountain  top  wind  speed 
Terrain  Roughness 

Global  Weather  Center  (GWC)  nomograms 


MOUNTAIN  WAVE  LLT 

Mountain  top  wind  speeds 
Cross  mountain  SLP  gradient 

Rate  of  decrease  of  Scorer  Parameter  with  height 
Harrison  nomogram 


LLT  INDICES 

burton's  Turbulence  Index  (3TI)  =  f(wind,  stability,  pressure 
tendency,  and  roughness) 

GWC  Mechanical  Turbulence  Index  =  f(wind,  roughness) 


Critical  values  for  the  various  parameters  listed  in  table  1  are  a  function  of 
the  geographical  area,  the  time  of  day  and  the  year,  and  the  aircraft 
category.  Many  values  are  listed  in  Lee  et  al,  (1979),  AWS  (1979),  and  FAA 
(1987).  The  most  frequently  quoted  lower  threshold  value  is  20  knots  for 
mountain  waves  and  low-level  mechanical  turbulence  significant  to  aircraft 
operation.  At  the  upper  end  of  the  scale,  50  knots  corresponds  with  severe 
turbulence  in  all  cases,  although  for  some  aircraft  the  threshold  is  sig¬ 
nificantly  lower  (35-40  knots).  Jones  et  al.  (1970)  have  completed  an  exten¬ 
sive  evaluation  of  the  use  of  wind  speed,  lapse  rate,  roughness,  BTI,  RI,  and 
Showalter  index  with  LO-LO  CAT  data.  Their  results  indicate  the  importance  of 
wind  speed,  roughness,  and  stability  and  the  utility  of  BTI  in  diagnosing  and 
predicting  LLT. 

For  the  most  part,  rules  of  thumb  for  LLT  forecasts/nowcasts  are  primarily 
numerical,  that  is,  related  to  the  parameters  listed  in  table  1.  Those  values 
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are  listed  in  the  previonsly  cited  references.  Several  rules  of  thumb  are 
available  for  specif''  localities;  they  are  usually  found  in  TFM.,  ale'ting 
local  forecasters  to  turbulence  prone  areas.  A  useful  summary  of  general 
rules  of  thumb  is  given  by  FAA  (ig??)  in  a  table  entitled  "Locations  of 
probable  turbulence  by  intensities  versus  weather  and  terrain  features."  Some 
other  useful  rules,  com, non  throughout  the  literature  in  one  form  or  another, 
rre  the  intensity  of  turbulence  always  increases  with  wind  speed  and  roughness 
and  the  turoulence  o‘leir>;nts  {hence  surface  gusts  and  LLT)  have  dimensions 
proportional  to  ruo  sire  of  the  roughness  elements. 

In  fiyurt}  I,  I oc.i  1 -tuni  ng  refers  to  the  procedure  of  adapting  the  cent''alized 
LLl  tore  .a .  1/ now.,  I t  to  Uie  local  area.  It  requires  a  careful  eval'i^..ian  of 
the  oentriiicoj  product.  Local  parameter  evaluation  (table  1),  PIKEP5,  and 
rules  0.  tiiumi;  are  introduced  to  tailor  the  forecast  to  the  needs  and  limita¬ 
tions  of  tie  lo'.a'!  u'-er.  Many  rules  of  thumb  are  specialized  for  a  particular 
1  oca le . 

"Me  t-wa  t'.ni  ng"  is  tne  co!nn'"i,  ter.n  for  monitoring  a  critical  situation  once  the 
forecds  c./nowcast  h..,s  been  made.  In  critical  evaluations  both  local -timing  and 
me  c-watchi  ng  ar;  labor-i  ntensi  va  (P.ichwien,  1979). 

Finally,  as  with  any  coi.iprehensi ve  forecast  scheme  a  systematic  verification 
is  carried  out  to  .ncni  tor  the  skill  and  improve  the  quality  of  LLT  forecasts 
(McGinley,  i9d6j;  this  step  is  often  unsatisfactory  because  of  the  few  verifi¬ 
cation  reports  from  aircraft. 

2.1,3  Future  Improvements  in  LLT  Forecasts/Nowcasts 

As  illustrated  in  tfia  jireceding  section,  current  procedures  of  LLT 
forecast!  ng/nowc  istiny  are  based  prirarily  on  the  establishi^d  relationship  of 
tiie  variou;  "types”  of  LLT  to  certain  large-scale  patterns  via  conceptual 
mesoscale  models  via  the  quantification  of  those  patterns  using  available 
data,  and  via  the  OKperience  and  attention  of  the  local  forecaster.  There  are 
several  obvious  probleii;^  with  this  scheme.  Although  the  synoptic  patterns 
associated  with  LLT  are  wel'  -known,  the  use  of  conceptual  mesoscale  models  is 
problematic.  Most  of  the  conceptual  models  currently  brought  to  bear  on  the 
operati.;nal  pr.)blen  are  "mean"  or  "typical"  two-dimensional  pictures  of 
phenomena  that  ivav.j  large  spatial  and  temporal  variabilities.  Therefore, 
tner.^  .vi  I  !  oe  nuny  situations  that  they  describe  poorly.  Also,  there  appears 
to  be  a  wide  variation  in  the  understanding  and  application  of  these  models  by 
forecasters.  Tno  fiattern  evaluation  step  can  overccme  some  of  these  problems, 
especially  wnon  tiia  t  step  is  accomplished  objectively.  Currently,  such 
ijuality  cuntr.i]  ,,in  by  assured  only  d*  a  weather  central.  Sooner  or  later  the 
1  Ll  forevst  will  reac.h  tne  local  forecaster  and  subjectivity  will  be  i  ntt  o- 
(juce  i  whe'i  tne  '.'■'jcast  is  tailored  to  the  local  area. 

Iho  final  Serious  or-jiilmn  that  plagues  LLT  forecasts/nowcasts  is  lack  of 
diti.  Lvn-n  ()l),ii"o I i v"  pattern  evaluation  procedures  suffer  from  the  lack  of  a 
very  '(.tenjiv'  ui'  uopi'M  s  ti  i,<i  ted  data  base.  Aside  from  the  use  of  a  few  semi- 
qu.i  nf  i  t !  l ;  Pii-ieFb.  none  of  the  current  procedures  is  based  on  direct  mea- 
b urein-:  n  tv  of  '  LT  . 

rnis  v:ct.ion  of  th  ■  report  describes  an  examination  of  the  literature  tfia  t  was 
perfor  ).‘d  fo  ].■  tn-’-'.'i  no  wimther  recent  research,  especially  in  the  areas  of  the 
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LLT  causes  discussed  above,  could  be  adapted  to  improve  LLT  forecasts.  Spe¬ 
cifically,  the  possibility  of  improving  the  model /parameter  components  of  the 
LLT  forecast  procedure  was  considered. 

During  the  last  10  years,  major  progress  has  been  made  in  the  areas  of  meso- 
scale  observations  and  mesoscale  modeling.  Thus  more  and  more  detailed 
descriptions  of  a  wide  variety  of  mesoscale  circulations  have  become 
avai lable. 

Two  recent  reviews  are  availaole  in  texts  on  mesoscale  circulations;  Pielke 
(1984),  which  covers  especially  modeling  studies  done  in  the  United  States  up 
to  the  early  ASCOT  papers,  and  Atkinson  (1981),  which  provides  a  somewhat 
broader  treatment  of  mesoscale  research  conducted  in  Europe  and  Asia  as 
well.  Both  texts  discuss  mathetfiatical  models  of  mesoscale  circulations; 
Atkinson  (1981)  also  summarizes  much  observational  data.  In  addition,  Ray 
(1986)  covers  much  the  same  material  from  a  forecasters  perspective.  Three 
overviews  of  "Mountain  Meteorology"  are  also  available:  GARP  (1978),  Smith 
(1979),  and  Heister  and  Pennel  (1980). 

Probably  the  most  usable  results  of  recent  mesoscale  research  are  those  from 
studies  of  moist  convection.  As  indicated  in  the  last  section,  much  of  the 
recent  inforination  on  gust  fronts,  outflow  boundaries,  downbursts,  and  related 
phenomena  have  been  adequately  reviewed  in  recent  works  by  Doswell  (1982, 
1985),  Fujita  (1985),  Ray  (1986),  and  many  others.  The  relative  ease  of 
observation  of  those  phenorivena  (compared,  for  example,  to  mechanical  turbul¬ 
ence)  and  their  role  in  several  fatal  and  we  1 1 -documented  aircraft  accidents 
(for  example,  Fujita,  1978,  1986)  likely  accounts  for  the  rapid  assimilation 
of  the  new  information  by  the  forecast  community.  This  is  not  the  case  for 
progress  made  in  other  areas  of  mesoscale  research. 

With  the  exception  of  the  use  of  satellite  imagery  to  locate  regions  of  moun¬ 
tain  lee  wave  activity  and  the  development  of  a  number  of  strong  downslope 
windstorm  (SOW)  prediction  aids  (for  example.  Brown,  1986),  few  dramatic 
improvements  have  been  made  in  the  prediction  of  mountain  waves  in  the  last  25 
years  (also  see  Uurran,  1986).  This  situation  exists  despite  an  intense 
research  effort  that  has  greatly  improved  our  understanding  of  those  phenomer... 
(for  example,  GARP,  1978;  Smith,  1979;  Heister  and  Pennel,  1980;  Mass  and 
Albright,  1985;  Kuettner,  1986). 

Information  that  has  evolved  from  research  but  still  awaits  application  to  the 
LLT  problem  includes  the  extension  of  the  simple  lee  wave  model  (Alaka,  1960), 
to  include  the  SOW  type  (Lilly  and  Zipser,  1972;  Lester  and  Fingerhut,  1974), 
a  better  understanding  of  the  dynamic  causes  of  SOW  (Klemp  and  Lilly,  1975, 
1978;  Peltier  and  Clark,  1979;  Smith,  1985;  Uurran,  1986d),  and  the  further 
use  of  satellite  data  to  diagnose  SOW  in  some  areas  (Elrod,  1986;  Lester  and 
Bach,  1986). 

Although  SOW  theory  has  been  advanced  significantly,  none  of  the  current 
models  have  been  adapted  to  the  prediction  of  the  details  of  SDWs  and  asso¬ 
ciated  LLT  in  an  operational  setting.  However,  there  are  some  parameters  from 
SOW  theory  that  may  be  useful  in  development  of  new  LLT  forecast  tools. 
Several  of  the  theoretical  studies  noted  above  have  emphasized  the  importance 
of  the  steepness  of  the  lee  slope  of  the  mountain  in  the  production  of  SDWs. 
Although  this  requirement  was  documented  many  years  ago  (for  example,  see 
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Harrison,  1966)  for  the  production  of  strong  lee  waves,  and  by  Scorer  (1978) 
and  many  others  for  the  separation  of  flow  over  a  ridge,  no  specific  applica¬ 
tion  was  encountered  in  the  LIT  prediction  methods  reviewed  to  date.  The 
ability  of  microcomputers  to  rranipulate  and  display  detailed  terrain  data 
bases  locally  indicates  that  the  time  is  right  for  the  local  forecaster  to 
make  good  use  of  that  information. 

Brinkmann  (1975)  and  Giusti  (1987)  nave  shown  that  the  stable  layer  which  is 
always  present  above  the  mountain  tops  during  mountain  wave  events  is  signif¬ 
icantly  stronger  during  SJWs.  Klemp  and  Lilly  (1975)  have  shown  that  if  the 
atmosphere  is  approximated  as  a  three-layer  hydrostatic  model,  and  disturbed 
by  a  'Dountain,  strong  surface  winds  will  be  produced  in  the  lee  in  proportion 
to  the  ampl  i  f  ication  factor 


AMP  =  Ni  X  N3/^N2)2 


where  Ni,  N2,  and  H?  are,  respectively,  the  Brunt-Vaisal la  frequencies  for  the 
stable  layer  at  and  below  mountain  top,  the  upper  troposphere,  and  the  strato¬ 
sphere.  Other  factors  such  as  tne  dimensions  of  the  layers,  the  vertical  wave 
lengths  of  the  lee  waves,  and  the  mountain  height  are  important,  but  current 
work  by  Lester,  Bach,  and  Muranaka  (1988),  suggests  that  AMP  evaluates  an 
important  contribution  of  SOW.  Giusti  and  MacKay  (1988)  are  currently  inves¬ 
tigating  the  predictive  value  of  AMP  via  regression  techniques. 

Recently,  Muranaka  (1988)  used  a  microcomputer  to  apply  the  iwss  consistent 
wind  model  developed  by  Ludwig  et  al.  (1985,  COHPLtX),  to  the  analysis  of  the 
surface  wind  distrioutions  during  two  Sows  over  the  foothills  of  the  Canadian 
Hockies.  The  results  showed  promise  for  operational  use,  and  it  was  recom¬ 
mended  that  further  experimentation  be  made  with  a  simple  two-dimensional  lee 
wave  model  as  the  upper  bounaary  to  determine  whether  the  COMPLEX  can  be  used 
pi  ac  tii-dl  ly ,  to  generate  nowcasts  of  surface  winds  and,  thus,  LLT. 

In  dealing  with  LLT  due  to  the  flow  of  stable  air  over  complex  terrain,  one  of 
tne  considerations  is  the  nature  of  eddies  downstream  of  individual  barriers. 
These  phenomena  include  hydraulic  jumps,  horizontal  meanders,  and  turbulence 
wakes  in  tne  lee  of  hills  and  other  terrain  obstruction.  The  Froude  number 
(f)  is  a  dimensionless  number  that  lies  in  certain  critical  ranges  when  the 
flow  takes  on  certain  unique  characteristics.  F  is  defined  as 


F  =  U(NL)-i  , 


where  d  is  the  undisturbed  fluid  velocity,  N  is  the  Brunt-Vai sal  la  frequency, 
and  L  is  a  character  i  sti  c  obstacle  dimension  (for  example,  the  height), 
liaines  (  1987)  nas  recently  given  an  extensive  review  of  the  interpretation  and 
application  of  the  Froude  number.  F  may  be  useful  for  LLT  diagnosis  in  com¬ 
plex  terrain  since  its  calculation  does  not  involve  any  assumptions  that  would 
be  invalidated  by  nonuni  form! ty  (unless,  of  course,  some  "bulk"  F  was  esti¬ 
mated  tor  a  large  area).  bincc  the  measure  is  taken  upwind  of  an  obstacle, 
local  instabilities  may  aliect  the  airflow  which  would  not  be  firedicted  by  the 
Froude  number.  Some  t/pical  studies  that  have  used  F  are  listed  in  table  II. 
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TABLE  2.  TYPICAL  STUDIES  THAT  HAVE  USED  THE  FROUDE  NUMBER 


Manins  and  Sawford  (1982)  Found  a  critical  Froude  number  of  1.6, 

above  which  the  synoptic  flow  would  flush 
the  small  valley  wliere  the  study  was  held. 
Concluded  that  tne  critical  Froude  number 
would  be  terrain  dependent. 


Furman  and  Wooldridge  (1984)  Calculated  Froude  number  for  flow  around 

and  over  an  obstacle.  Found  a  value  of  F  = 
0.09  for  very  stable  flow  around  the  hill. 
For  flow  that  just  begins  to  go  over  the 
hill  an  F  =  0.4  was  calculated;  for  smooth 
flow  over  the  hill  F  >  2.0.  The  authors 
conclude  that  the  flow  over  the  hill  at  the 
lower  Froude  number  was  caused  in  part  by  an 
unstable  region  of  tne  windward  base  of  the 
hill  which  did  not  affect  the  Froude  number. 


.(Ooldridge  and  Furman  (1984)  Observations  of  a  simple  hill  and  flow 

parameterized  oy  Froude  number.  For  values 
0.3  <  Fr  <  0.7  superpressured  balloons 
passed  around  the  hill  and  were  occasionally 
caught  in  a  lee-side  rotor  which  was  not 
present  at  higher  Froude  numbers. 


bmi th  (1984)  Calculated  Froude  numbers  for  observed 

flows.  Conclusions  tentative. 


As  discussed  earlier,  the  treatment  of  nechanical  turbulence  in  current  LIT 
forecasts  does  not  depend  on  a  conceptual  model  beyond  a  simplistic  idea  of 
random  eddies  that  become  stronger  as  the  wind  increases  and/or  the  terrain 
becomes  rougher.  Ihe  improvement  of  instrumentation  coupled  with  the  increase 
in  the  study  of  the  transport  and  diffusion  of  pollutants  nas  lead  to  a  much 
better  understanding  of  the  small-scale  flow  features  that  develop  near  land- 
sea  boundaries,  and  complex  terrain.  i4any  details  are'  found  in  Atkinson 
(1981)  and,  wi  tli  respect  to  numerical  modeling,  in  PieUe  (1984).  The  spe¬ 
cific  details  for  individual  studies  vary  widely  as  a  t unction  oi  the  ter¬ 
rain.  examples  of  some  of  those  details  are  shown  i  ri  table  3. 
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TABLE  3.  MECHANICAL  TURBULENCE  STUDIES  DETAILS 


I'^hrer  and  Pielke  (1978) 
dornstein  and  Tnoripson  (iydl) 
Yosnikado  (idol) 

Lyons  (197k;) 

Lyons  (1975) 

McCarthy  and  Young  (1978) 

i-lanis  (t  Sawford  (  1979) 
Dickerson  &  Cudiksen  (1981) 
Netf  and  King  (  1985) 

Porch  et  al .  ( 1988) 

5 tone  &i  Hoard  (  1988) 

Whiteman  (1988) 

Whiteman  (  1982) 

Whiteman  &  i'lcKee  (  1982) 


Sea  breeze  with  onshore  synoptic  winds. 

Sea  breeze  front  retardation. 

Synoptic  scale  influence  on  sea  breeze. 

Climatology  and  prediction  of  sea  breeze. 

Turbulent  diffusion  at  shoreline. 

Profile  (structure)  of  New  England  coastal 
front.  Authors  note  that  the  front,  once 
understood,  is  quite  predictable.  This 
could  well  be  true  of  many  local  phenomena. 

Model  of  katabatic  winds. 

ASCOT  report  (Geysers  studies). 

Studies  using  acoustic  sensors  (ASCOT). 

Contributions  of  valley  tributaries. 

Side  wall  circulations  and  flow  surges. 

Vertical  profiles  of  downslope. 

Observed  vertical  and  cross  valley 
structure  using  tethered  balloons. 

Observed  inversion  breakup  and  provided 
time  estimates  based  on  inversion  structure. 


Segal  et  al.  (1986) 

Wooldridge  <14  Orgill  (19/8) 

Selvain  et  al.  (1983) 

Baker  et  al .  (1984) 

Etling  and  Wamser  (1988) 


Observations  and  modeling  of  the  effect  of 
cloud  shading  on  sea-breeze  and  upslope 
wi nds. 

Momentum  flux  over  mountain  valley,  observa¬ 
tions  of  synoptic  scale  flow  penetrating  the 
top  of  a  va  1 1  ey  . 

Discussion  of  the  mechanism  by  which  rotor 
circulations  are  maintained. 

Importance  of  angle  of  ridgeline  to  flow. 

Structure  of  vortexes  in  lee  of  islands. 
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Since  one  of  the  causes  of  LIT  is  dry  convection,  it  nas  .jeen  recommended  by 
some  scientists  that  the  application  of  some  of  the  basic  boundary  layer 
stability  criteria  such  as  Pasquill  staoility  class,  kichardson  Number,  and 
Moni n-Obukhov  Length  (for  example,  Hunn,  1966;  Haugen,  1973)  be  considered  as 
LLT  forecast  parameters.  Pasquill  noted  that  tne  downwind  spread  of  a  plume 
or  puff  was  dependent  on  the  distance  from  the  source  ind  the  stability  of  the 
atmosphere  in  the  area  as  long  as  one  had  some  knowledge  of  the  behavior  of 
the  atmosphere  in  general.  In  order  to  treat  such  boundary  layer  influences 
objectively,  he  (Pasquill  1961)  developed  a  set  of  diffusion  curves.  The 
curves  were  derived  from  experiments  carried  out  over  relatively  smooth  roll¬ 
ing  terrain.  Stability  conditions  are  estimated  from  "surface"  windspeed  and 
daytime  insolation  (a  very  general  measure  of  the  temperature  profile).  For 
each  stability  class  a  curve  is  drawn  relating  tne  horizontal  or  vertical 
dispersion  of  a  scalar  quantity  to  the  distance  from  the  source. 

txtension  of  this  kind  of  study  into  complex  terrain  was  complicated  by  tne 
fact  that  mountain  valleys  tended  to  form  extremely  stable  layers  during  the 
night  and  that  the  daytime  mixing  generally  exceeded  ttiat  estimated  by  the 
Pasquill  curves  for  similar  conditions  in  flat  terrain  (for  example,  Koch 
etal.,  1977).  In  addition,  small-scale  convergence  zones  and  return  flow 
situations  made  the  estimation  of  diffusion  much  mor-  difficult.  Even  with 
the  application  of  quantitative  means  to  estimate  staijility,  this  categoriza¬ 
tion  of  turbulence  has  not  proven  to  be  valid  for  complex  terrain  where  there 
is  often  a  gradient  in  stability  and  wind  velocity  and  therefore  in  turbul¬ 
ence. 

It  should  also  be  noted  that  for  air  pollution  concerns,  once  the  atmosphere 
is  "well  mixed"  there  is  no  further  need  to  categorize  turbulence.  Since, 
according  to  the  Pasquill  stability  classes,  this  occurs  at  a  windspeed  of 
6  meters/second,  in  our  opinion  these  categories  have  little  bearing  on  tur¬ 
bulence  encountered  during  aircraft  flight,  even  for  light  aircraft. 

The  flux  Richardson  Number  (R^)  is  defined  as 


where  is  the  eddy  diffusivity  for  heat,  is  the  eddy  diffusivity  momen¬ 
tum,  and  R^  is  the  gradient  Richardson  number  given  by 


R.  =  9-ii 
’  0  dz 


In  the  last  equation,  g  is  gravity,  Q  is  potential  temperature,  and  u  is 
windspeed.  Richardson  (1920)  derived  tne  expression  for  R^-  for  "just 
turbulent"  flow  in  horizontally  homogeneous  conditions.  may  bo  related  to 
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tne  Moni n-Obukhov  scale  length,  L  (Monin  and  Obukhov,  1953;  Panosky  and 
Dutton,  1984)  as 


z/P.  ,  R.  >  0 
L  = 

(l-5R.)z/R.,  Ri  >  0 


The  flow  is  better  characterized  by  L  since  R^  is  a  function  of  height  (z). 
The  similarity  theory  is  that  the  boundary  layer  expands  or  contracts  with  L 
(Munn,  I960).  During  the  daytime,  the  ratio  (-z/L)  represents  the  relative 
importance  of  heat  convection  {-z/L  strongly  negative),  mechanical  turbulence 
(z/L  approximately  zero),  and,  at  night,  the  suppression  of  mechanical  turbul¬ 
ence  (z/L  strongly  positive).  Thus  the  function  z/L  yields  more  information 
aoout  the  type  of  turbulence  over  its  entire  range  of  values.  R^  on  the  other 
hand  simply  exceeds  some  finite  critical  value  for  the  onset  of  the  turbul¬ 
ence.  The  application  of  either  L  or  R^-  to  the  characterization  of  LLT  on  an 
operational  basis  is  questionable  because  of  the  poor  quality  of  observations 
normally  available.  Furthermore,  in  areas  of  complex  terrain,  the  representa¬ 
tiveness  of  these  variables  is  in  doubt. 

Although  many  studies  and  models  of  mesoscale  circulations  exist,  few  are 
operationally  useful  predictive  numerical  models.  Most  are  diagnostic,  and 
while  they  can  increase  understanding  of  the  mechanisms  that  drive  the  cir¬ 
culations  they  do  not  provide  real-time,  operational  forecasts.  Prognostic 
models  such  as  the  work  done  by  Yamada  (1981,  1983)  are  very  detailed  and,  in 
order  to  make  reasonable  predictions,  would  require  an  excessive  amount  of 
computer  time.  As  several  authors  have  concluded,  there  is  no  numerical  model 
for  complex  terrain  today  that  is  a  good  forecasting  tool.  Table  4  shows  the 
works  consulted  to  reach  this  opinion. 

One  developing  use  of  numerical  models  is  the  application  of  simple  mass 
consistent  wind  models  to  diagnose  wind  distributions  in  complex  terrain  on 
the  basis  of  a  few  observations.  Although  the  validity  of  such  models  still 
awaits  testing,  their  requirement  for  only  small  computer  resources  as  well  as 
some  promising  results  from  preliminary  experiments  (see  next  section)  sug¬ 
gests  that  their  application  should  be  pursued. 

2.1.4  The  Application  of  AI/Expert  Systems  to  the  LLT  Forecast  Problem 

AI  is  a  generic  term  referring  to  the  use  of  a  computer  to  imitate  hunan 
oenavior  that  is  generally  thought  to  require  intelligence.  Expert  systems 
(LS)  and  knowledge-based  systems  (KBS)  are  less  stringent  terms  dealing  with 
the  use  of  computers  to  emulate  human  thought  processes  under  stricter  guide¬ 
lines  (using  empirical  relationships  based  on  experience  and  knowledge  of  the 
programmer).  Racer  and  Laffney  (1984)  introduced  the  term  interpretive 
processing  (IP)  as  an  application  of  ES/KBS  in  meteorological  applications  and 
quote  tne  following  definitions. 
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TABLE  4.  EXAMPLES  OF  MESOSCALE  NUMERICAL  MODELS 


Kosoerg  ct  al.  (1576) 
Dickerson  (1978) 

Sherman  (1978) 

Erasmus  (1986a) 

Yarnada  (1981) 

Yamada  (1983) 

Meyers  et  al .  ( 1985) 

Yamada  and  Kao  (1986) 

Me  11  or  and  Yamada 

Wyngaard  (1985) 

Wyngaard  (1985d) 

Erasmus  (1986b) 

Henmi  (1988) 


A  simple  mass  consistent  ivind  model. 

Mass  consistent  wind  model  specifically  for 
mountainous  terrain. 

Mass  consistent  wind  model. 

Multilayer  mass  consistent  model  of  Oahu  under 
Trade  Wind  influence. 

Complex,  9-level,  model  of  nocturnal  drainage 
flows  including  a  multilevel  soil  moisture  and 
heat  flux  model 

Description  of  a  simplified  turbulence  model 
that  is  still  highly  complex  relative  to  a 
simple  diagnostic  model. 

A  Large  Eddy  Simulation  (LES)  model.  (This  is 
essentially  a  work-in-progress  report.) 

Simulation  of  the  hiarine  boundary  layer  during 
CATE,  3-D,  3nd  moment,  turbulence  closure 
model . 

(1974)  Comparison  of  models  for  turbulence  in  the 
planetary  boundary  layer. 

Considers  in  general  terms  the  value  of  current 
modeling  efforts  and  suggests  that  there  may  be 
room  for  cost-effective  improvement. 

Describes  in  general  terms  LES  techniques  and 
suggests  where  this  form  of  predictive  numerical 
model  might  prove  useful. 

Evaluation  of  a  diagnostic  inass-consistent 
model  against  observed  data. 

Compares  a  complex  multilayer  model  with  a 
simplified  model.  Found  that  the  simplified 
model  showed  unrealistic  windspeeds. 
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AI  is  a  subfield  of  computer  science  concerned  with  the  concepts  and  methods 
of  symbolic  inference  by  a  computer  and  the  symbolic  representation  of  the 
knowledge  to  be  used  in  imking  inferences.  A  computer  can  be  made  to  behave 
in  ways  that  humans  recognize  as  "intelligent"  behavior  in  each  other 
{Feigenbaum  and  McCorduck,  1983). 

Al  is  the  development  of  computer  programs  that  can  solve  problems  normally 
thought  to  require  huirvan  intelligence  (Uuda  and  Shortliffe,  1983). 

ESs .. .[are] .. .problem-sol vi ng  computer  programs  that  can  reach  a  level  of 
performance  comparable  to  that  of  a  human  expert  in  some  specialized  problem 
domain  (Nau,  1983). 

...a  K8S  is  an  AI  program  whose  performance  depends  more  on  the  explicit 
presence  of  a  large  body  of  knowledge  than  on  the  possession  of  ingenious 
computational  procedures,  by  expert  system  we  mean  a  KBS  whose  performance  is 
intended  to  rival  that  of  human  experts  {Duda  and  Shortliffe,  1983). 

IP  is  defined  as  a  computer  interactive  procedure  that  enhances  the  abilities 
of  the  weatner  forecaster  to  decide  on  a  forecast.  The  procedure  makes  it 
easier  to  draw  conclusions  from  the  meteorological  analysis  of  observational 
data,  forecasting  techniques,  and  past  forecaster  experience  available  when 
deciding  on  a  forecast. 

The  possible  applications  of  AI  to  meteorology  cover  a  spectrum,  ranging  from 
decision  trees,  such  as  developed  -by  Brown  (1986)  to  forecast  programs  capable 
of  learning  (baffney  and  dacer,  1983)  and  beyond.  The  National  Weather 
Service  (NWS)  is  increasing  its  automation  of  field  operations  as  part  of  its 
modernization  efforts,  with  one  of  its  areas  of  concentration  being  the  field 
of  IP.  Since  the  forecast  problem  involves  reduction  of  available  data, 
identification  of  significant  data  and  guidance  (numerical  and  manual),  and 
the  application  of  both  explicit  and  implicit  relationships,  rules  of  thumb, 
etc.  to  create  a  forecast  product,  a  competent  IP  system  would  be  of  great 
benefit.  Racer  and  Gaffney  (1984)  give  an  example  of  a  prototype  IP,  further¬ 
more,  they  (Racer  and  Gaffney,  1984)  envision  a  three-fold  benefit  from  the 
application  of  K8S/ES  to  weather  forecasting 

(1)  to  provide  improved  data  analysis  and  decision-making  support  due  to 
enhanced  consistency  and  thoroughness, 

(2)  to  support  training  of  new  forecasters, 

(3)  to  support  skill  maintenance  for  experienced  forecasters,  especially 
with  regard  to  their  actions  in  infrequently  occurring/unfamiliar  situations. 

Unquestionably,  an  LET  forecast  system  that  accomplished  the  above  items  would 
go  far  in  alleviating  the  LET  forecast  problem  at  regional  and  local  forecast 
offices  (for  exam()le,  see  next  section). 

AI  technology  is  being  used  in  varying  degrees  as  a  forecast  tool.  As  noted 
Brown  (1986)  has  developed  a  simple  decision  tree  approach  to  forecasting 
downslope  windstorms  in  Colorado.  His  is  a  program  using  "if-then"  structures 
to  consolidate  significant  data  (both  analysis  and  numerical  guidance)  into  a 
valid  indicator  of  the  probability  of  strong  downslope  winds.  Gaffney  and 
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Kacer  (1983)  have  developed  a  prototype  system  for  severe  storm  advisories 
that  is  capable  of  "learning"  behavior.  This  system  is  based  upon  formalized 
rules  developed  by  Crisp  (1959)  and  Miller  (1972)  of  the  Air  Force  GWC.  Racer 
and  Gaffney  (1986)  quote  a  personal  communication  with  J.  T.  Schaefer  of  the 
National  Severe  Storms  Forecast  Center  detailing  a  KBS  that  includes  "a  severe 
weather  checklist  of  10  parameters  that  are  evaluated  as  a  group  using  "if 
then"  rules  to  determine  the  “possibility"  of  a  storm."  Racer  and  Gaffney 
(1986)  also  detail  a  diagnosis  procedure  for  evaluating  numerical  guidance 
materials  developed  by  Simpson  (1971)  at  the  NWS  National  Hurrican  Center.  It 
used  a  decision  ladder  to  systematically  analyze  the  performance  of  numerical 
models  with  the  goal  of  improving  them. 

There  is  an  apparent-gap  in  the  spectrum  of  technical  applications  of  AI  to 
weather  forecasting.  Gaffney  and  Racer's  "learning"  program  is  at  the  high 
end,  but  it  is  only  a  prototype.  The  checklist/decision  tree  approach  (for 
example,  Lee,  1988),  at  the  low  end  of  the  spectrum,  is  the  only  application 
of  AI  commonly  in  use.  While  this  is  an  improvement  over  manual  methods,  much 
greater  benefits  could  be  realized  by  the  use  of  "smarter"  systems. 

Possible  candidates  would  be  the  refinement  of  the  dependence  of  a  local 
turbulence  index  (LTI)  (see  next  section)  on  wind  direction,  relative  weight¬ 
ing  of  input  parameters,  etc.  The  system  would  likely  use  numerical  input/ 
output  rather  than  a  "natural  language"  user  interface  characteristic  of  full 
AI  systems.  The  program  should  include  a  training  mode  with  blocks  for; 
(1)  general  forecaster  familiarization,  (2)  training  on  use  of  program,  and 
(3)  practice  cases/case  studies. 

LLT  forecasting/nowcasting  would  greatly  benefit  from  the  application  of  AI 
concepts,  especially  with  the  greater  availability  of  powerful  microcomputers 
and  reasonably  priced  remote  sensing  devices.  In  an  effort  to  solve  the  LLT 
prediction  problem  at  Fort  Irwin,  Lee  (1988),  has  recently  developed  and 
tested  an  interactive  LLT  forecast  aid  to  be  operated  in  parallel  with  a  local 
observation  network  and  an  objective  wind  analysis  program  (Henmi  et  al., 
1988).  Results  are  promising.  In  the  current  study  (next  section),  an  LLT 
forecast  method  based  on  a  modification  of  one  of  the  indices  presented  in 
table  1  is  proposed  for  use  with  a  local  data  base  and  an  efficient  objective 
analysis  technique  to  determine  the  local  wind  field.  As  will  be  seen,  the 
proposed  technique  could  be  easily  adapted  for  forecaster  interaction  and 
continuous  nowcast  display,  two  important  attributes  of  a  practical  AI  system. 

In  addition  to  the  literature  review  above,  a  limited  search  was  conducted  for 
data  bases  that  include  simultaneous  tower  and  aircraft  measurements  of  LLT. 
Such  data  could  allow  the  development  of  improved  statistical  forecast  tech¬ 
niques  of  LLT.  By  far,  the  most  comprehensive  LLT/tower  data  base  collected 
to  date  was  LO-LO  CAT  (Jones  et  al.,  1970).  Data  from  Project  PHOENIX,  a  more 
recent  data-collecting  program,  offer  promise.  It  is  currently  being  analyzed 
by  Lilly  et  al .  (1988).  National  Center  for  Atmospheric  Research  (NCAR) 
maintains  many  data  bases  from  its  Research  Aviation  Facility  (RAF)  aircraft 
investigations,  however,  those  have  not  been  examined  closely.  Personal 
communications  with  interested  scientists  generally  discouraged  the  search  and 
the  attempted  utilization  of  such  data  bases  as  expensive,  with  a  good 
probability  of  being  unproductive.  Further  investigation  is  needed. 
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It  is  clear  from  tne  previous  sections  that  a  dramatic  improvement  in  LLT 
forecasti  ng/nowcdsti  ng  awaits  the  wide  availability  of  sophisticated  boundary 
layer  measurements  as  well  as  realistic  mesoscale  numerical  models  that  can  bo 
operated  economically  in  near  real  time.  Since  these  improvements  are  not 
likely  to  oe  seen  by  the  operational  forecaster  for  many  years,  a  practical 
recommendation  is  that  current  methods  be  used  more  efficiently.  On  the  basis 
of  the  literature  review  and  the  Fort  Irwin  problem  discussed  below,  it  is 
clear  that  significant  improvements  in  local  LLT  forecasts/nowcasts  can  be 
made  by  combining  the  simplest  AI  approaches  (for  example,  checklist/decision 
tree)  with  available  forecast/nowcast  techniques  (table  1)  and  with  current 
computer,  communication,  and  measurement  technology. 

2.2  The  LLT  Problem  at  Fort  Irwin 

2.2.1  Background 

LLT  associated  wi  tn  high  winds  at  the  NTC  at  Fort  Irwin,  California,  occa¬ 
sionally  causes  aircraft  accidents  or,  more  often,  the  cancellation  of  mis¬ 
sions  of  helicopters  operating  in  support  of  training  exercises.  The  latter 
problem  occurs  because  the  helicopters  may  not  fly  into  areas  where  severe 
turbulence  is  observed  or  predicted.  High  winds  and  turbulence  are  common 
during  fall  and  spring  and  proportionately  less  missions  should  be  expected  to 
be  flown  during  those  periods.  However,  there  is  a  general  perception  (no 
statistics  available)  among  command  and  flying  personnel  at  NTC  that  incorrect 
"overforecasts"  of  severe  LLT  contribute  significantly  to  the  total  number  of 
cancel  led  missions. 

It  is  the  purpose  of  this  portion  of  the  current  study  to  document  the  Fort 
Irwin  LLT  problem  and  to  seek  answers  to  the  questions:  "Can  the  LLT  fore¬ 
casts  oe  improved  significantly?  and  if  so,  how?" 

2.2.2  Procedures 

The  Fort  Irwin  proolem  was  documented  by  means  of;  (1)  a  review  of  the  NTC 
Forecaster  Handbook  (1987),  (2)  a  tour  of  the  Post,  (3)  interviews  with  the 
Science  Advi sor ,  local  pilots,  and  forecasters,  (4)  a  review  of  George  Air 
Force  Base  forecasting  procedures,  (5)  an  examination  of  a  number  of  past  case 
studies  of  LLT  incidents  on  or  near  the  Post,  and  (6)  a  questionnaire  com¬ 
pleted  by  NIC  permanent  party  pilots. 

2.2.3  Some  Important  Characteristics  of  LLT  at  Fort  Irwin 

Fort  Irwin  occupies  an  area  of  about  30  x  30  miles^  in  the  Mojave  Desert  35 
miles  northeast  of  Barstow,  just  south  of  Death  Valley.  It  lies  about  85 
miles  to  the  north  of  the  San  Bernardino  and  San  Gabriel  Mountains  and  about 
tne  same  distance  to  the  east  of  the  southern  Sierra  Nevada  where  the  highest 
peaks  exceed  10,000  feet  mean  sea  level  (m.s.l.).  On  the  Post,  the  terrain  is 
characterized  uy  rugged  peaks  separated  by  broad  valleys  (figure  2).  tleva- 
tions  (m.s.l.)  over  the  Post  vary  widely  from  near  6,100  feet  in  the  northeast 
corner  of  i,300  (feet)  in  the  southeast. 
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Ft.  Irwin  Topography  (Heights  in  Feet)  MSL 
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Figure  2.  Fort  Irwin  topography. 

LIT  that  causes  the  most  difficult  mission  scheduling  problems  at  Fort  Irwin 
during  fall  and  spring  is  caused  mainly  by  interactions  of  the  rugged  terrain 
with  very  strong  winds  that  occur  during  the  passage  of  fronts,  and  with  moun¬ 
tain  waves  (NTC  Forecaster  Handbook,  1987).  A1  thougli  significant,  convec- 
tively  produced  LIT  occurs  throughout  the  year,  it  does  not  present  a  problem 
as  serious  as  the  frontal  and  mountain-wave  generated  LIT.  This  is  evidently 
due  to  a  better  forecaster/pilot  understand! ng  of  convective  phenomena. 
Furthermore,  many  of  the  most  important  characteristics  of  moist  convection 
are  easily  identifiable  by  eye  and/or  by  radar.  Also,  there  is  a  strong 
relationship  between  thermal  activity,  time  of  day,  and  specific  terrain 
features.  In  contrast,  the  production  of  mechanical  turbulence  is  not  as  well 
understood  by  pilot  or  forecaster;  it  occurs  on  smaller  time  and  space  scales 
and  is  further  complicated  by  the  extreme  variations  in  terrain  across  Fort 
Irwin.  For  these  reasons,  the  study  of  the  Fort  Irwin  LIT  problem  concen¬ 
trates  on  mechanically  produced  LIT. 

2.2.4  Turbulence  Cases 

Synoptic  conditions  were  examined  for  cases  where  aircraft  encountered  signif¬ 
icant  turbulence  on  or  near  Fort  Irwin.  Dates  of  the  13  incidents  (furnished 
by  the  U.S.  Army  Atmospheric  Sciences  Laboratory  (ASL))  are  listed  in  table 
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5.  Although  other  cases  were  available  (Lee,  1988),  they  are  not  considered 
here  because  either  the  dates  fell  outside  the  October-May  period  (two  cases) 
or  there  was  insufficient  inforniation  (three  cases).  The  analysis  of  tne 
remaining  cases  included  the  documentation  of  synoptic  patterns  (location  of 
fronts,  troughs,  examination  of  wind,  contour,  pressure,  and  pressure  tendency 
gradients,  etc.),  unfortunately,  more  detailed  analysis  was  only  possible  for 
two  cases  because  the  turbulence  incidents  were  poorly  documented,  that  is, 
the  location  and  ti'ne  and,  in  some  i nstances,  the  date  could  not  be  verified. 
These  uncertainties,  coupled  with  tne  already  sparse  distribution  of  surface 
and  radiosonde  stations  in  the  area,  prohibited  further  systematic  analysis. 


TABLE  5.  DATES  OF  SIGKIFICAMT  LET  INCIDENTS* 


ia  Apr  7t) 
14  Apr  83 
lo  Jan  84 
7  tlov  84 
27  I’Tar  85 
7  Oct  85 
0  I'lar  8o 


♦furnished  by  ASL 


25  Nov  86 
6  Jan  87 
16  Jan  87 
5  Feb  87 
5  Mar  87 
19  Apr  87 


The  primary  result  or  trie  case  study  macroanalysis  was  that  9  of  tiie  13  cases 
occurred  in  the  vicinity  of  surface  fronts;  one  was  associated  with  a  sharp 
upper  level  trough  (no  clear  surface  front);  and  one  was  related  to  a  surface 
high-pressure  system  located  over  the  Great  Basin.  These  conditions  agree 
with  what  is  known  generally  about  LLT  turbulence-producing  processes  at  Fort 
Irwin  (NTC  Forcaster  Handbook,  1987);  strong  surface  winds  (and  therefore  LLT) 
are  caused  either  by  mountain  waves  or  by  strong  pressure  gradients  in  the 
vicinity  of  cold  fronts.  Furthermore,  these  features  are  identifiable  in  the 
large-scale  synoptic  pattern  and  should  therefore  be  anticipated  by  trained 
forecasters  using  standard  analysi  s/"met-watch"  procedures. 

2.2.5  Pilot  (juestionnaires 


As  noted  in  tne  literature  review  and  verified  above,  synoptic  patterns  are 
good  indicators  of  the  probability  of  LLT  occurrence  over  a  broad  area. 
However,  the  documentation  of  details  of  the  distribution  of  LLT  over  an  area 
of  the  size  of  Fort  Irwin  clearly  requires  subsynoptic  scale  information.  The 
fulfillment  of  tnis  requirement  is  difficult  because  historical  weather  data 
(pressure,  temperature,  winds,  etc.)  for  the  Post  are  primarily  limited  to  a 
single  location,  the  Bicycle  Lake  Army  Airfield  (8YS).  Furthermore,  aside 
from  the  poorly  documented  cases  listed  in  table  5,  data  on  turbulence  occur¬ 
rences  around  the  Post  are  nonexistent. 

In  order  to  overcome  this  problem,  at  least  partially,  a  questionnaire  related 
to  LLT  was  developed  and  distributed  to  Fort  Irwin  pernement  party  helicopter 
pilots.  That  group  was  selected  (for  example,  as  opposed  to  rotation  pilots) 
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because  of  their  interest  in  the  problem,  knowledge  of  the  Post,  and  similar 
and  extensive  flying  experience.  The  14  pilots  who  responded  to  the  question¬ 
naire  had  an  average  accumulated  flying  time  of  1946  hours,  with  an  average  of 
847  hours  at  Fort  Irwin. 

The  questionnaire  instructed  the  pilots  to  describe  LIT  not  associated  with 
thunderstorms.  They  were  also  asked  to  characterize  typical  missions,  to 
document  sources  of  meteorological  information,  to  evaluate  forecasts,  and  to 
give  suggestions  for  forecast  improvement.  A  sample  questionnaire  is  given  in 
appendix  A;  a  summary  of  all  resporises  except  those  related  to  the  spatial 
distribution  of  LIT  are  contained  in  appendix  b.  The  responses  to  the  LLT 
distribution  questions  are  described  briefly  below  and  presented  in  detail  in 
appendix  C. 

Figure  3  summarizes  the  major  result  of  the  questionnaire,  that  is,  an  esti¬ 
mate  of  the  location  of  the  primary  LLT  problem  areas  on  the  Post.  It  is  a 
composite  of  individual  maps  prepared  oy  the  pilots  in  response  to  the 
instruction:  "on  the  attached  map,  circle  those  locations  which  have  the 

highest  frequency  of  turbulence  significant  to  your  operations."  The  associa¬ 
tion  of  those  areas  with  particular  topographical  features  is  noted.  Many  of 
the  pilots  who  responded  to  the  questionnaire  provided  detailed  comments  about 
specific  turbulent  areas  (numbered  areas  in  figure  3).  The  comments  are  given 
in  appendix  C.  These  should  prove  useful  in  tne  future  for  both  pilot  and 
forecaster  training.  Altnough  the  utility  of  such  infonnation  is  clear, 
caution  is  advised.  The  sample  is  small  ahd  the  pilot's  reports  are  biased 
towards  the  primary  operational  areas;  thus  the  map  is  simply  a  guide  and  does 
not  represent  a  "climatology"  of  any  sort. 

2.2.6  Causes  of  the  LLT  Forecast  Problem  at  NTC 

Un  the  basis  of  the  site  visit  and  the  questionnaires  described  above,  it  is 
clear  that  the  primary  cause  of  the  LLT  problem  at  Fort  Irwin  originates  with 
the  production  of  a  point  forecast  for  NTC  by  the  weather  detachment  at  George 
Air  Force  Base  (which  carries  the  primary  responsibility  tor  local  LLT  fore¬ 
casts  for  Fort  Irwin).  Despite  the  fact  that  part  of  the  Fost  will  often  have 
severe  turbulence  conditions  as  predicted.  Fort  Irwin  personnel  will  perceive 
such  forecasts  as  erroneous  "overforecasts"  because  a  gradient  in  turbulence 
intensity  exists  across  tne  Post,  that  is,  a  significant  portion  of  NTC  is 
still  "flyable"  although  the  entire  Post  is  closed. 

The  problem,  therefore,  is  not  always  one  of  inaccurate  forecasts,  but  rather 
that  the  user's  needs  have  exceeded  the  forecaster's  current  capabilities. 

The  perception  of  “bad"  forecasts  has  damaged  forecaster  credibility  and  has 
led  to  forecast  "shopping"  by  pilots  and  the  necessity  for  command  personnel 
to  override  forecasts  in  order  to  accomplish  missions. 

Another  unfortunate  ramification  of  the  perceived  problem  is  that  rotation 
forecasters  at  Bicycle  Lake  (many  of  whom  are  not  familiar  with  NTC  forecast¬ 
ing  problems)  do  not  benefit  greatly  from  their  NTC  experience.  They  are 
often  left  out  of  the  forecast  loop,  relegated  to  observer  or  briefer  status, 
and  ignored  by  both  pilots  and  command  personnel.  Based  on  our  perception, 
they  are  essentially  an  under-utilized  resource. 
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Although  the  forecast  methods  used  Oy  VCV  to  produce  lLT  predictions  for  Fort 
Irwin  reflect  currently  accepted  tnetiiods  used  oy  the  Air  Weather  Service  (for 
example,  AwSP  idb-bb,  19y6)  and  other  agencies,  tlieir  gross  detail  cannot 
solve  tile  Fort  Irwin  problem.  In  order  to  solve  tne  problem  described  above, 
a  new  forecast  system  is  needed,  that  is,  one  that  recognizes  the  scale  of  the 
turbulence  and  the  required  detail  of  the  forecasts. 

A  minimal  but  practical  LLT  forecast/ nowcast  system  would  include  (1)  a  meso- 
scale  network  of  surface  observation  stations  across  the  Post  with  data  acces¬ 
sible  in  real  time,  (2)  a  system  of  regular  pilot  reports  of  turbulence 
intensity  over  the  Post,  and  (3)  an  objective  analysis  (nowcasting)  and  fore¬ 
cast  system  based  on  uoth  large  and  local  scale  information.  The  surface 
observation  stations,  coupled  with  local  LLT  reports  would  provide  the  data 
base  for  the  development  and  continued  improvement  of  objective  nowcast/ 
forecast  techniques. 

In  addition  to  the  three  items  listed  above,  there  is  an  important  fourth 

component  for  a  successful  system... (4)  the  human  factor.  As  was  clearly 
noted  in  the  literature  review,  successful  forecasting  over  small  space  and 
time  scales  requires  "met-watchi ng,"  that  is,  it  is  a  labor-intensive  task. 
The  rotation  forecaster  at  Fort  Irwin,  given  the  proper  techniques  and  train¬ 
ing,  can  certainly  bear  a  significant  part  of  the  LLT  forecast  responsibil¬ 
ity.  The  only  drawback  foreseen  will  be  maintaining  continuity  from  rotation 
to  rotation. 

Assuming  that  items  (1),  (2),  and  (4),  above,  will  be  available,  tne  next 

section  discusses  the  development  of  a  potential  method  for  the  objective 
analysis  (nowcasting)  and  forecasting  system  for  Fort  Irwin. 

1.2,1  An  LLT  Index-Based  Forecast  Scheme 

A  useful  objective  forecast  scheme  for  LLT  must  evaluate  the  contributions  of 
all  significant  physical  processes  that  cause  the  turbulence.  It  must  be 
practical,  that  is,  easy  to  learn,  easy  to  use,  and  based  on  available  data 

bases.  As  discussed  in  the  literature  review,  the  3TI  nas  been  used  as  a 

major  input  to  worldwide  LLT  predictions  by  GWC  of  the  LI.S.  Air  Force  during 
the  mid-sixties  and  early  seventies  (Burton,  1964;  Burnett,  19/0).  It  was  a 
successful  objective  technique  (for  example,  Jones  et  al.,  1970),  generally 
satisfying  the  requirements  stated  above  and  specifically  addressing  terrain 
roughness.  In  the  following,  BTI  is  described  in  detail,  a  modification  is 
proposed  to  adapt  it  to  an  area  the  size  of  Fort  Irwin,  and,  finally,  an 
example  of  its  use  is  demonstrated  for  a  case  of  severe  LLT. 

Ihe  BTI  (for  example,  Burnett,  1970)  is  given  by 


BTI=iT  +  V  +  5  +  T,  (1) 


where  K  is  "roughness,"  the  difference  between  the  highest  and  lowest  terrain 
features  in  the  area  of  interest  in  hundreds  of  feet.  V  is  windspeed,  at  2000 
feet  above  ground  level  (AGL)  in  knots.  S  is  "stability,"  the  lapse  rate 
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(10  X  Utb  C/iOUU  teot)  in  the  lowest  100  millibars.  T  is  "synoptic  forcing" 
absolute  value  of  tno  J-hour  pressure  tendency  in  tenths  of  a  millibar. 

Threshold  values  of  ill  for  various  turbulence  intensities  and  category  1 
aircraft  are  given  in  taole  6. 


TABLE  6. 

BTI  AS  A  FUNCTION  OF 

TURBULENCE 

LEVEL 

Light 

bU 

Modera  te-Severe 

90 

Severe 

100 

Modera  te 

70 

Extreme 

120 

ttTI  has  been  successfully  applied  to  large  areas  (for  example,  the  Mojave 
Desert),  however,  an  aoaptation  to  an  area  of  the  size  of  NIC  causes  several 
difficulties.  Altnougti  the  values  of  wind,  stability,  and  pressure  tendency 
may  be  similar  for  tne  Si.ialler  area,  the  numerical  value  of  the  roughness 
component  will  usually  be  less  than  its  value  for  a  larger  area.  It  follows 
that  total  bTI  values  v.'iil  also  be  smaller  with  roughness  being  weighted  less 
for  the  smaller  area.  Thus,  threshold  BTI  values  for  critical  LLT  will  be 
lower  tnan  tiiose  definea  by  past  experience  (table  2).  Another  problem  that 
arises  is  that  observations  of  winds,  stability,  and  pressure  tendency 
Observations  are  not  usually  available  on  the  same  scale  as  the  terrain 
information.  In  orde'^  to  deal  with  these  problems,  the  BTI  computation  has 
been  modified  in  the  tallowing  manner.  First,  it  is  assumed  that  the 
stability  (b)  and  ^)r(?saur.■  tendency  T  primarily  reflect  large-scale  influ¬ 
ences,  and  nay  oe  reer^senie J  uy  rather  gross  measurements,  for  example,  a 
single  sounding  and  on.-  rojCf-sonta  ti  ve  pressure  tendency  for  all  scales  of 
importance.  Locally,  iiu^ever,  roughness  (R)  and  windspeed  (V)  may  differ 
widely  from  their  irauni.'i,  ]■_  Ijes. 

Tne  modified  index  v.i/np,>  ta  ti  jn  is  then  performed  in  two  steps.  Initially  a 
nacroscale  BTI  (Lq.  (i)i  is  defined  for  the  broad  area  encompassing  Fort 
Irwin.  That  index  is  Lien  '.o.iitied  for  local  measurements  of  roughness  and 
wind  to  yield  a  "Local  scale  Turbulence  ^ndex,"  that  is. 


lTI.  =  oTI  X  C.R.  +  Vj/[R  +  V), 

1  II  ^ 


(2) 


where  BTI  is  given  by  eg.  (i),  suscript  i  represents  one  of  several  measure¬ 
ment  points  or  gridpoints,  and  subscript  x  indicates  the  maximum  measured 
gridpoint  value.  In  order  to  maximize  LTI,  in  the  case  that  the  denominator 
of  Lq.  (2)  exceeds  oil,  then 

LTIi  =  [R  +  V).  .  (3) 
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2.2.8  The  Application  of  LTl 

The  LTI  computation  was  tested  to  determine  whether  clear  gradients  of  the 
index  would  develop  across  the  Post  under  large-scale  conditions,  which  are 
diagnosed  as  turbulence  as  indicated  by  BTI.  Since  no  good  data  base  is 
available,  computations  were  made  with  data  from  a  known  case  of  significant 
LLT.  The  case  chosen  was  tfie  27  March  1985  incident  involving  the  crash  of  an 
0-2  at  Fort  Irwin  during  mountain  wave/LLT  conditions. 

Since  the  distribution  of  winds  was  not  available,  surface  (iU  in)  winds  were 
interpolated  at  2  kilometer  grid  points  by  means  of  COMPLEX  adapted  to  the 
Fort  Irwin  terrain.  Input  consisted  of  terrain  data  at  the  grid  points  and  a 
simple  wind  profile  based  on  surface  and  850  millibar  data  estimated  for  BYS. 

For  the  computation  of  BTI,  macroscale  roughness  was  determined  as  a  maximum 
terrain  height  difference  across  the  Post.  For  LTI,  "local"  roughness  at  grid 
points  was  estimated  from  topographical  charts  as  the  maximum  difference  in 
terrain  height  surrounding  each  grid  point.  The  3-hour ly  pressure  tendency 
was  determined  as  an  interpolated  value  from  surrounding  weather  stations, 
large-scale  stability  was  computed  from  the  estimated  surface  to  850  millibar 
lapse  rate  at  BYS,  and  the  macroscale  for  BTI  was  estimated  from  850  millibar 
at  BYS. 

Since  the  purpose  of  the  experiment  was  to  determine  whether  clear  LTI  gra¬ 
dients  would  develop  under  realistic  inacroscale  conditions,  the  boundary  layer 
structure  of  COMPLEX  was  arbitrarily  set  to  produce  a  maximum  response  in 
surface  winds  with  the  given  wind  profile.  Calculations  of  LTI  were  iTiade  for 
the  time  of  the  turbulence  incident  (1500Z)  and  6  hours  before  and  after. 

The  LTI  analyses  are  presented  in  figures  4,  5,  and  6  for  eacl)  time  period.* 
Despite  the  high  values  of  BTI  (also  see  table  2).  LTI  varies  widely  across 
the  Post  in  each  case.  As  would  be  expected,  there  is  a  strong  spatial  corre¬ 
lation  with  terrain  features  (compare  with  figure  2)  indicating  a  strong 
contribution  by  local  roughness.  If  these  are  realistic,  they  suggest  that 
LTI  may  offer  help  to  the  forecaster  for  the  discrimination  between  "flyable" 
and  "unflyable"  areas  across  the  Post  on  occasions  when  the  macroscale  fore¬ 
cast  would  close  tne  entire  area  to  flying. 

The  inspection  of  LTI  analyses  in  figures  4,  5,  and  b  emphasizes  other  poten¬ 
tial  advantages  of  the  index.  For  example,  it  is  easy  to  interpret  objec¬ 
tively  by  both  forecasters  and  pilots,  lending  itself  to  a  simple  display  (for 
example,  on  a  monitor).  COMPLEX  and  similar  objective  wind  analysis  schemes 
run  quickly  on  available  microcomputers  and  thus  such  analyses  and  the  index 
could  be  updated  in  nearly  real  time  and  displayed  continuously.  Furthermore, 
predicted  values  of  BTI  from  large-scale  forecasts  could  be  updated  locally 
with  local  values  of  V  and  K, 


*Analyses  of  COMPLEX  windspeed  and  wind  direction  for  each  time  are  presented 
i n  appendi x  U, 
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Figure  4.  LTI  for  0900Z. 
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Figure  5.  LTI  for  1500Z. 
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The  major  shortcoming  of  the  LTI  computation  presented  here  is  that  detailed 
local  wind  observations  are  nonexistent.  The  validity  of  the  scheme  described 
here  depends  strongly  on  the  availability  of  wind  observations  and  on  the 
accuracy  of  the  wind  interpolation  scheme.  It  should  be  pointed  out  that  the 
computation  of  low-level  turbulence  index  (LLTI)  at  the  site  of  the  wind 
observations  is  not  dependent  on  the  interpolation  algorithm.  The  success  of 
COMPLEX  interpolations,  on  the  other  hand,  is  strongly  dependent  on  knowledge 
of  the  shape  of  the  top  of  the  boundary  layer.  Furthermore,  COMPLEX  does  not 
explicitly  deal  with  nonhydrostatic  phenomena  that  are  often  characteristic  of 
strong  wind  regimes  in  rugged  terrain.  Several  other  questions  await  a 
suitable  turbulence  and  surface  data  base.  These  include  the  following: 
(1)  what  is  the  smallest  horizontal  scale  on  which  LTI  can  be  calculated  and 
still  be  meaningful?  {2  kilometers  were  selected  arbitrarily  in  the  present 
case);  (2)  what  is  the  relative  importance  of  each  variable  (R,  S,  T,  V)  as 
the  scale  decreases?  (3)  what  is  the  possibility  of  modifying  LLTI  as  a 
function  of  wind  direction  to  take  into  account  lee  waves  and  other  wake 
phenomena?  (4)  would  an  interpolated  wind  at  a  higher  level  be  more  appro¬ 
priate  for  LTI  computation  than  a  10-meter  wind? 

2.2.9  Development  of  an  Expanded  Data  Base 

The  previous  sections  have  demonstrated  that  significant  improvement  of  LLT 
forecasts  at  Fort  Irwin  will  require  a  data  base  consisting  of  surface  obser¬ 
vations  and  turbulence  observations  on  the  scale  of  the  desired  forecasts. 
The  literature  review  has  revealed  that  data  bases  of  the  desired  quality  are 
rare,  in  general,  and  nonexistent  for  Fort  Irwin,  specifically.  Current  plans 
to  place  a  numoer  of  automated  remote  weather  stations  around  the  Post, 
coupled  with  the  success  with  the  pilot  questionnaire  discussed  earlier, 
suggest  that  the  establishment  of  a  program  to  regularly  acquire  and  archive 
pilot  report  (PIREPS)  could  provide  such  a  data  base.  In  the  short  run,  such 
a  program  would  open  up  forecaster/pilot  communication  channels  on  a  more 
regular  basis.  In  the  long  run,  the  data  base  would  provide  information  for 
the  systematic  development  of  improved  LLT  nowcast/forecast  methods  on  a 
small-scale  desired  by  pilots  and  command  personnel. 

Towards  this  end,  an  efficient  pilot  report  form  was  developed  and  distributed 
to  Fort  Irwin  pilots.  A  sample  of  the  form  is  presented  in  appendix  E.  Data 
are  currently  being  archived  at  ASL,  White  Sands  Missile  Range,  New  Mexico. 


3.  SUMMARY  OF  RESULTS,  CONCLUSIONS,  AND  RECOMMENDATIONS 

The  literature  review  has  revealed  that  current  LLT  forecast  techniques  are  a 
result  of  the  need  to  forecast  a  small-scale  phenomenon  with  large-scale  data. 
Thus  aside  from  occasional  PIREPs,  the  forecast  techniques  lean  heavily  on 
large-scale  pattern  recognition  and  evaluation  of  those  patterns  with 
available  large-scale  data.  Although  the  current  forecast  techniques  lend 
themselves  to  automation  and  objectivity  at  the  weather  central,  at  the  local 
level,  a  large  amount  of  subjectivity  is  introduced  as  a  function  of  the  local 
forecasters  experience.  Furthermore,  by  the  very  nature  of  LLT,  its  accurate 
prediction  is  labor  intensive  at  the  local  level,  and  the  manpower  is  not 
always  present,  especially  in  critical  situations. 
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The  literature  also  indicated  that  improvements  in  LLT  in  the  near  future 
would  most  likely  come  from  the  improvement  of  local  data  acquisition  and 
objective  analysis  and  display,  that  is,  utilizing  current  instrument 
communications  and  microcomputer  technology. 

The  LLT  problem  at  Fort  Irwin  is  a  good  example  of  the  problem  with  current 
LLT  forecasting.  Although  standard  LLT  forecast  methods  are  used  and  there  is 
no  evidence  that  the  resulting  accuracy  is  not  what  should  be  expected,  the 
perception  of  the  forecast  users  is  that  the  forecasts  are  inaccurate.  The 
reason  for  this  perception,  is  that  the  users'  needs  simply  exceed  the  current 
forecast  capabilities. 

A  survey  of  the  turbulence  problem  was  completed  primarily  as  a  forecast  aid. 
A  prototype  objective  forecast  scheme  was  then  developed  to  eliminate 
variations  in  forecast  quality  while  capturing  the  smaller  scale  details  of 
the  turbulence  field  as  required  by  NTC  operations.  Finally,  a  local  PIREP 
form  was  developed  to  increase  forecaster  knowledge  of  the  problem  and  to 
provide  a  data  base  for  the  development  of  better  forecasts. 

The  literature  review  was  not  comprehensive.  For  example,  further  investiga¬ 
tion  of  the  results  of  Atmospheric  Studies  in  Complex  Terrain  (ASCOT)  and  the 
Alpine  experiment  (ALPtX)  should  be  conducted  to  determine  applications  to 
LLT.  Also,  there  remains  a  serious  need  to  develop  techniques  to  diagnose  and 
forecast  the  presence  of  significant  low-level  wind  shear  in  elevated  inver¬ 
sions,  especially  in  complex  terrain.  Further  work  also  needs  to  be  done  on 
the  development  and  validation  of  BTI  for  LLT.  The  current  data  gathering  at 
Fort  Irwin  offers  a  unique  data  base  for  this  purpose,  and  its  analysis  should 
be  pursued. 
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APPENUIX  A.  PILOT  QUESTIONNAIRE 


Pilot  Questionnaire 

1.  Pilot  Experience 
Name 

Rank 

Career  flying  time  (hours) 

Flying  time  at  Ft.  Irwin  (hours) 

Flying  time  in  current  aircraft  type(s) 

2.  Helicopter  Operating  Limitations 

List  current  aircraft  type(s)  and  their  operating  limitation 
with  respect  to  turbulence,  surface  wind  conditions,  and 
visibility. 

3.  Typical  Mission 

If  possible,  characterize  a  "typical  mission:” 

a)  Base  of  operations 

b)  Most  common  take  off  time(s) 

c)  Are  there  any  times  of  day  or  night  when  you  never  fly? 

d)  Most  of  your  flights  occur  below  what  altitude  (indicate 
ASL  or  AGL) ? 

e)  What  is  the  length  of  most  of  your  missions  (hours  and 
tenths) ? 

f)  In  what  area(s)  of  Ft.  Irwin  are  most  of  your  missions 
flown? 
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4.  Characteristics  of  Turbulence  Mot  Associated  With  Thunderstorms 

We  are  concerned  with  Low  Level  Turbulence  MOT  associated  with 
thunderstorms.  With  that  in  mind,  please  answer  the  following 
questions : 

Fill  out  the  table  below  with  respect  to  Ft.  Irwin 


Low  Level  Turbulence 

Time  of  Day 

Season 

General  Weather  Pattern 

Worst  Problem 

Least  Problem 

5.  Low  Level  Turbulence  Areas 


On  the  attached  map,  circle  those  locations  which  have  the  highest 
frequency  of  turbulence  significant  to  your  operations.  Clearly 
number  each  location  you  have  circle  on  the  map  and  list  that 
number  below,  giving  a  brief  description  of  the  turbulence  problem 
including  the  wind  direction  when  the  turbulence  problem  occurs . 


6.  Turbulence  Incidents 

List  below  the  dates,  times,  and  locations  of  any  particularly 
notable  turbulence  incidents  or  accidents  at  Ft.  Irwin.  If  neces¬ 
sary,  circle  the  location  on  the  attached  map  and  label  clearly 
with  a  number  fro  cross  referencing.  Give  as  many  details  as 
possible. 

7.  LLT  Forecasting  Accuracy 

Comment  on  the  following  specific  aspects  of  low  level  turbulence 
forecasts  for  Ft.  Irwin. 


Accuracy  of  turbulence  intensity  forecasts  (adequate,  over-,  or 
under-forecast? ) 


Accuracy  of  forecasts  of  location,  areal  extent  of  turbulence 
(over-,  or  under-forecast?) 
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8.  Preflight  Briefings 

Where  do  you  get  your  preflight  briefings?  When  do  you  get  them 
(relative  to  take  off  time)? 

9.  Inflight  Advisories 

Do  you  ever  receive  inflight  weather  (turbulence)  advisories?  If 
yes,  from  where? 

Do  you  ever  give  PIREPS?  If  yes,  to  whom? 

10.  Postflight  Debriefings 

Do  you  ever  give  post  flight  weather  debriefings?  If  yes ,  to 
whom? 

11.  Recommendations  to  improve  LLT  forecasts 

What  are  your  recommendations  to  improve  low  level  turbulence 
forecasts  for  your  operations  at  Ft.  Irwin? 
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APPENDIX  B.  SUMMARY  OF  PILOT  QUESTIONNAIRES 


Summary  of  Pilot  Questionnaires 

In  order  to  better  define  the  LLT  problem  at  Ft.  Irwin,  a  quei>tionnaire 
was  distributed  to  the  permanent  party  helicopter  pilots  at  the  site 
visit  to  Ft.  Irwin  in  October,  1987.  Fourteen  questionnaires  were 
returned  by  Ft.  Irwin  Pilots  and  are  summarized  below. 

1.  Pilot  experience  (13  responses) 


Mean 

Median 

Range 

Career 

flying  time  (hrs) 

1946 

1500 

800-5400 

Flying 

time  at  Ft.  Irwin  (hrs) 

847 

700 

300-1880 

Flying 

time  in  type(s) 

UH-IH  (12  responses) 

1720 

700 

500-5400 

Other  insufficient  responses 

2.  Helicopter  operating  limitations 

Surface  wind  conditions  ("to  crank")  : 

Maximum  Hind  30  knots  (35  knots  for  OH-58?) 

Maximum  Gust  Spread  15  knots 
Turbulence : 

May  not  fly  into  region  of  reported  or  predicted  severe 
turbulence  (exception:  waiver  by  Post  Commander) 

Flight  not  recommended  into  areas  where  moderate 

turbulence  has  been  reported  by  Category  2  or  higher 
aircraft. 

Ceiling/Visibility:  Restrictions  in  uncontrolled  mountainous 

terrain.  Day:  0.5  mi/500  ft;  night:  1.0  mi/1000  ft. 
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3.  "Typical  Mission"  (#  responses) 

a)  Base  of  Operations:  Barstow/Dagget  (9);  Field  Site,  BYS 
(4) 

b)  Most  conmon  take-off  tines:  0300-0700L  (13) 

c)  Times  of  day  or  night  when  flight  is  restricted  None  (8); 
Might  (5) 

d)  Most  flights  conducted  below:  1000  ft  (1)  ;  500  ft  (1)  ; 

300  ft  (7);  200  ft  (3);  100  ft  (2) 

e)  Mission  duration:  Mean  (9):  3.1  hrs;  Duty  day  (6): 

12  +  hrs 


f)  Area(s)  of  most  missions:  Entire  Post  (9);  Central  and 
Southern  Corridors  (3);  South  (1);  BYS  (1) 


Characteristics 

of  LLT  and  not  associated  with 

thunderstorms . 

Time  of  Day 

Season 

Weather  Pattern 

Worst 

problem 

Noon-evening  (12) 

Sumner  (3) 

Near  front  (2) 

Morning  (4) 

Fall/Spring 

(8) 

Strong  (sw) 

winds  (5) 

Winter  (1) 
all  (2) 

Other  (3) 

Least 

problem 

Morning  (7) 

Winter  (4) 

Light  winds  (3) 

Night  (2) 

Fall/Spring 

(2) 

Other  (3) 

Other  (2) 

Summer  (2) 
All  (3) 

5.  See  Appendix  C 

6 .  See  Appendix  C 

7.  Comments  on  Accuracy  of  LLT  forecasts  for  Ft.  Irwin 
Intensity 

Fcsts  of  LLT  intensity  too  strong  (6)  ;  fcsts  ok  but  tending  to  be 
too  strong  (2)  ;  fcsts  ok  but  tending  to  be  too  weak  (1)  .  Three 
respondents  simply  indicated  that  LLT  intensity  fcsts  were  inaccu¬ 
rate. 

Area 

Fcsts  of  LLT  area  too  large  (9);  fcst  area  ok  (1);  area  too  small 
(2) ;  three  respondents  simply  indicated  that  fcst  areas  were 
inaccurate . 


8.  Preflight  briefings  are  obtained  within  1-1.5  hrs  of  takeoff  from 
one  of  the  following  sources:  DAG  FSS ,  VC7,  BYS. 
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9. 


Inflight  advisories  are  obtained  from  BYS  (8);  DAG  FSS  (2);  other 
acft  (1);  never  (4).  Two  responses  indicated  that  BYS  is  not 
accessible  from  the  air  for  such  information.  One  respondent 
stated  that  BYS  dees  not  pass  on  inflight  advisories. 

All  14  pilots  indicated  that  they  gave  PIREPS  at  one  time  or 
another  to  BYS  (12),  to  other  acft  (1),  and  to  DAG  FSS  (3). 
Although  not  stated,  the  conflict  between  these  numbers  and  those 
in  the  paragraph  above  indicate  that  some  of  these  PIREPS  are  not 
given  in  flight,  but  are  given  after  the  acft  have  landed. 

10.  Postflight  weather  debriefing  have  been  given  at  various  times  by 
nine  (9)  of  the  respondents  to  BYS  (4)  ,  Unit  Operations  (4)  ,  DAG 
FSS  (2)  ,  VCV  (1)  .  Five  (5)  respondents  never  give  post  flight  de¬ 
briefings  . 

11.  Recommendations  to  improve  low  level  LLT  forecasts  at  Ft.  Irwin: 

a)  Permanent  party  forecaster  at  Ft.  Irwin  (10) 

b)  Install  remote  surface  wind  sensors  (9) 

c)  More  PIREPS  (2) 

d)  Other:  Better  communications,  turbulence  recovery  route  to 
Barstow/Daggett ,  permanent  party  observer 
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APPENDIX  C.  RESPONSES  TO  PILOT  QUESTIONNAIRE 


Responses  to  Pilot  Questionnaire  on  Low  Level  Turbulence 
Distribution  at  Ft.  Irwin 
(Questionnaire  Itens  #5  and  #6) 

5 •  Low  Level  Turbulence  Areas  (non-convective) 

On  the  attached  map,  circle  those  locations  which  have  the  highest 
frequency  of  turbulence  significant  to  your  operations.  Clearly 
number  each  location  you  have  circled  on  the  map  and  list  that 
number  below  giving  a  brief  description  of  the  turbulence  problem, 
including  wind  direction  when  the  turbulence  problem  occurs. 


6 .  Turbulence  Incidents  (non-convective) 

List  below  the  dates,  times,  and  locations  of  any  particularly 
notable  turbulence  incidents  or  accidents  at  Ft.  Irwin.  Give  as 
many  details  as  possible. 


RESULTS 


5.  The  map  on  the  following  page  is  a  composite  of  those  loca¬ 
tions  which  have  the  highest  frequency  of  turbulence  significant 
to  helicopter  operations  at  Ft.  Irwin  (based  on  11  responses) . 

The  Roman  Numerals  on  the  map  correspond  with  descriptions  of  the 
turbulence  in  the  respective  areas  described  in  Table  1.  The  most 
frequently  listed  areas  were  Tiefort  Mountain  (8)  ,  Granite  Pass 
(9)  ,  and  the  area  just  East  of  Drinkwater  Lake  (7)  . 

6.  None  of  the  respondents  gave  complete  descriptions  (e.g.,  in¬ 
cluding  dates,  times,  and  locations)  of  any  particularly  notable 
turbulence  incidents  or  accidents  at  Ft.  Irwin.  Some  partial 
descriptions  are  listed  in  Table  Cl. 
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TABLS  Cl 


I  GRANITE  MTNS 

II  BUNKER  - 


III  AVAHATZ  MTNS 


IV  BNASITB  - 


1.  Hind  from  5H  10  knots  or  greater 

2.  Strong  up  and  down  drafts 

3 .  Ridge  wave 

1.  Hind  H  to  SH  below  200  ft.  AGL 

2.  East  of  Drinkwater  bake 

3.  Hinds  from  Hest  about  40  knots 

4.  Air  flow  thru  valley  (Westerly)  can  cause 
strong  winds  and  turbulence 

5.  Hinds  W  to  NW  severe 

1.  Pinacle  approaches  wind  -  South 

2.  Hind  changes  direction  up  to  30°  wind  out 
of  S/SW  Oct.  1982  during  downward  portion 
of  an  approach  to  a  6000’  pinnacle  the  UH- 
1  was  caught  in  a  severe  downdraft,  then 
updraft  causing  shoulder  harness  to  lock 
eat.  wind  velocity  20-30  knots 

1.  Hind  from  Hest,  wind  shifts  30° 
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TABLE  Cl 
(Continued) 


V 

TIEFORT  MTNS  - 

1. 

Sudden  tailwind  gusts.  Wind  up  to  100° 
difference  between  landing  zones. 

2. 

MTN  wave 

3. 

Strong  updrafts  and  downdrafts  -  occurs 
when  the  wind  picks  up 

VI 

COYOTE  CANYON  - 

1. 

Heavy  turbulence  when  elsewhere  cal.Ti. 

Winds  from  West 

2. 

Has  moderate  turbulence  at  times  of 
weather  warnings 

VII 

SOUTHWALL  - 

1. 

Hinds  tend  to  come  from  the  rear  of 
hovering  helicopters.  Downdraft 

VIII 

CHINAMAN'S  HAT  - 

1. 

Swirling  winds 

2. 

Wind  out  of  S/SW  Sept.  1986,  AUH  -1,  Fully 
loaded  encountered  light  to  moderate 
turbulence  while  in  slow  flight  through 
saddle.  The  aircraft  was  in  a  20-30  knot 
tailwind  condition  -  causing  power 
applications  resulting  in  momentary  loss 
of  rotor  RPA. 

IX 

TRAINING  AREA  B  - 

1. 

Winds  tend  to  be  stronger  here  than  other 
parts  of  the  Post. 

X 

GRANITE  PASS  - 

1. 

Turbulence,  wind  gusts 

2. 

Pass  turbulence 

3. 

Moderate  turbulence  at  time  of  weather 
warnings 

4. 

Saddle  causes  venturi  effect 

5. 

Winds  W  to  NW  severe 

6. 

Moderate  to  severe  short  duration  surface 
winds  from  the  North 

XI 

RED  PASS  - 

1. 

Turbulence  and  wind  gusts  make  flying 
difficult 

XII 

LEACH  LAKE  PASS  - 

1. 

Tailwind,  gusts,  downdrafts 

2. 

Winds  out  of  SW  Feb  -  Mar  87  during 
straight  and  level  flight  along  ridge 
lines,  the  OH55  aircraft  encountered 
moderate  turbulences.  Est  wind  velocity 

XII 

SOUTHWALL  - 

1. 

Wind  tends  to  come  from  the  rear  of 
hovering  helicopters  downdrafts 

XIV 

1. 

During  sustained  high  wind  conditions,  is 
subjecc  to  turbulence  due  to  the  Venturi 
effect. 

XV 

EAST  GATE  - 

1. 

Winds  W  to  NW  severe 

XVI 

CHECK  POINT  - 

1. 

Winds  W  to  NW  severe 

XVII 

HIDDEN  VALLEY  - 

1. 

Winds  W  CO  NW  severe 

XVIII 

BIKE  LAKE  - 

1. 

Winds  W  to  NW  severe 

XIX 

1. 

Ridge  Winds 

XX 

1. 

During  sustained  high  wind  conditions,  is 
subject  to  turbulence 
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APPENDIX  D.  COMPLEX  MINOSPEEDS  AND  DIRECTIONS  FOR  21  MAR  85 


COMPLEX  Nindspeeds  and  Directions  for 
21  Mar  85  at  0900Z,  1500Z  and  2100Z 
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APPENDIX  E.  POST-FLIGHT  TURBULENCE  SURVEY 


DATE. 


PQST-FUGHT  TURBULENCE  SURVEY 


EXAMPLE 


ON  THE  MAP  ON  THE  SEVEHSE.  REPORT  THE  FOLLOWING;  OR 

AIRCRAFT  ID  NOT  REQUIRED'  BE  Ao  SPECIFIC  AS  POSSIBLE. 

1.  LOCATION. 

2.  TII-IE  (LOCAL). 

.3.  ALTITUDE  ABOVE  GROUND  LEVEL  (AGL)  . 

4.  TURBULENCE  INTENSITY.  USE  STANDARD  REF-DETING  CATEGORIES. 
None(  N) ,  Lijsht(L/ , ModerateCM :  ,3ev»re{S)  .eXtr^mefX/ 

PLEASE  INCLUDE  NEGATIVE  (N)  REPORTS. 

T.  SEE  EXAJ1PLE  ABOVE. 


BIBLIOGRAPHY 


Abbs,  D.J.  and  R.A.  Pielke,  1986:  Thermally  forced  surface  flow  and 
convergence  patterns  over  northern  Colorado.  Mon.  Weath.  Rev., 
114,  2281-2296. 

Alaka,  M.  A.,  1960:  The  airflow  over  mountains.  WMO  TN  34,  Geneva, 
135  pp. 

Atkinson,  B.W.,  1981:  Meso-Scale  Atmospheric  Circulations.  Academic 
Press,  London,  495  pp. 

AWS ,  1979:  Meteorological  Techniques,  (with  Revisions).  AWSP  105-56. 

AWS ,  1985:  Catalog  of  Air  Weather  Service  Technical  Documents. 
AWS/TC-85/001,  133  pp. 

Bains,  P.,  1987:  Upstream  Blocking  and  the  Air  Flow  over  Mountains, 
Annual  Review  of  Fluid  Mechanics ,  19 ,  75-97 . 

Baker,  R.W.  et  al . ,  1984:  Ridgecrest  winds  in  mountainous  terrain. 

Presented  at  the  Third  Cqnferencj*  on  Mountain  Meteorology,  Ameri¬ 
can  Meteorological  Society,  Boston,  Mass. 

Becker,  R.J.  and  K.B.  Gage,  1972:  Limited  shear  zone:  a  hazard  to 

aviation.  Presented  at  the  International  Conference  on  Aerospace 
and  Aeronautical  Meteorology,  American  Meteorological  Society, 
Boston,  Mass. 

Bedard,  A.J.,  1982:  Source  and  detection  of  atmospheric  wind  shear. 

A I AA  Jour.  ,  2_p  ,  940-945. 

Biggs,  W.G.  and  M.E.  Graves,  1962:  A  lake  breeze  index.  J.  Appl. 

Meteor . ,  1,  474-480. 

Bonner,  W.D.  and  J.  Paegle,  1970:  Diurnal  variations  in  boundary 

layer  winds  over  the  south  central  United  States  in  summer.  Mon 
Weath._Rev_.  ,  98,  735-744. 

Bornstein,  R.D.  and  W.T.  Thompson,  1981:  Effects  of  frictionally 

retarded  sea  breeze  and  synoptic  frontal  passages  on  sulfur  diox¬ 
ide  concentrations  in  New  York  City.  J.  Appl.  Meteor.,  20,  843- 
858. 

Bowen,  B.M.,  J.M.  Dewart,  and  A. I.  Chen,  1983:  Stability  class 

determination;  a  comparison  for  one  site.  Presented  at  the  Sixth 
Symposium  on  Turbulence  and  Dif fusion,  American  Meteorological 
Society,  Boston,  Mass. 

Bradbury,  T. ,  and  J.  Kuettner,  (Ed),  1976;  Forecasters  Manual  for 
Soaring  Flight,  OSTIV  Scientific  Section,  Geneva,  119  pp. 


57 


Brinkmann,  W. ,  1975:  Strong  downslope  winds  at  Boulder,  Colorado. 

MWRr  102,  592-602. 

Brown,  J.  M. ,  1986,  A  Decision  Tree  for  Forecasting  Severe  Downslope 

Windstorms  in  Colorado,  Eleventh  Conference  on  Weather  Forecasting 
and  Analysis,  (AMS),  17-30  Jun  86. 

Buettner,  K.J.K.  and  N.  Thyer,  1966:  Valley  winds  in  the  Mt.  Rainier 
area.  Arch.  met.  geophys.bigklim.j_B,  14,  125-147. 

Burnett,  P.,  1970:  Turbulence  Forecasting  Procedures,  USAF  AFGWC  TM 
70-7.  78  pp. 

Burton,  M. ,  1964:  Unpublished  Notes. 

Businger,  J.  A.,  1973:  Turbulent  transfer  in  the  atmospheric  surface 
layer.  Atmos.  Environ.,  14,  741-750 

Camp,  D.,  and  W.  Frost,  1987:  Atmospheric  Turbulence  relative  to 

aviation  missile  and  space  programs.  NASA  Conference  Pub  2468, 

265  pp. 

Caracena,  F.,  R.A.  Maddox,  L.R.  Hoxit,  and  C.F.  Chappell,  1979: 

Mesoanalysis  of  the  Big  Thompson  storm.  Mon.  Weath.  Rev,,  107,  1- 
17. 

Colson,  D.,  1954:  Meteorological  problems  in  forecasting  .mountain 
waves.  Bull .  _Amer._  Meteor,.  __Soc.,  35,  363-371. 

Colton,  D.E.,  1976:  Numerical  simulation  of  the  orographically  induced 
precipitation  distribution  for  use  in  hydrologic  analysis.  J. 

Appl. _ Meteor.,  15,  1241-1251. 

Coulman,  C.E,,  J.R.  Colquhoun,  R.K.  Smith,  and  K.  Mclnnes,  1985: 

Orographically  forced  cold  fronts  -  mean  structure  and  motion. 
Boundary-Layer  Met . ,  32,  57-83. 

Crisp,  C.  A.,  1959,  Training  Guide  for  Severe  Weather  Forecasters, 
Document  AFGWC/TN-79/COL,  Air  Weather  Service  (MAC),  Air  Force 
Global  Weather  Central,  Nov  59. 

Defant,  F.,  1950:  Theorie  der  land-  und  seewind.  Arch  Meteorol. 

Gegphys.  Bioklimatol.  Ser  A,  2,  404-425. 

Defant,  F.,  1951:  Local  Winds.  "Compendium  of  Meteorology."  American 
Meteorological  Society,  Boston,  Mass.,  pp.  655-672. 

Dickerson,  M.H,,  1978:  MASCON-A  mass-consistent  atmospheric  flux  model 
for  regions  with  complex  terrain.  J.Appl.  Meteor.,  17,  241-253. 

Dickerson,  M.H.  and  P.H.  Gudiksen,  1981:  "ASCOT  FY-1981  Progress 

Report",  Lawrence  Livermore  National  Laboratory  Report  UCID-13878- 
81,  ASCOT-2. 


58 


Doswell,  C.,  1982:  The  operational  meteorology  of  convective  weather. 
Volume  I:  Operational  mesoanalysis ,  NOAA  TM  NWS  NSSFC-5,  95  pp . 

Doswell,  C.,  1985:  Volume  II:  Storm  scale  analysis,  NOAA  TM  ERL  ESG- 
15,  240  pp. 

Duda,  R.  O. ,  and  E.  H.  Shortliffe,  1983:  Expert  Systems  Research, 
Science,  Vol.  220,  15  Apr  83,  p  261. 

Durran,  D. ,  1936a:  "Mountain  Waves"  in  Mesoscale  Meteorology  and 
FqrecasJ:ing ,  (P.  Ray,  ed.)  AMS,  472-492. 

Durran,  D. ,  1986b:  Another  Look  at  Downslope  Windstorms.  Part  I:  The 
development  of  analogs  to  supercritical  flow  in  an  infinitely  deep 
continuously  stratified  fluid.  JAS ,  43,  2528-2543. 

Ellrod,  G. ,  1986:  Uses  of  Satellite  Imagery  in  Turbulence  Detection. 

Paper  AIAA-86-Q258 ,  AIAA  24th  Aerospace  Sciences  Meeting,  January 
6-9,  Reno.  5  pp. 

Erasmus,  D.A,,  1986(a):  A  model  for  objective  simulation  of  boundary- 
layer  winds  in  an  area  of  complex  terrain.  J.  Climate  Appl . 
Mjetepr,,  25,  1832-1841. 

Erasmus,  D.A.,  1986(b):  A  comparison  of  simulated  and  observed 

boundary-layer  winds  in  an  area  of  complex  terrain.  J.  Climate 
Appl .  Meteor .  ,  25,  1842-1852. 

Estoque,  M.A. ,  1962:  The  sea  breeze  as  a  function  of  prevailing 
synoptic  situation.  J.  Atmos.  Sci.,  19,  244-250. 

Etling,  D.  and  M.  Wamser,  1988:  Some  aspects  of  atmospheric  vortex 
streets  formed  in  the  wake  of  large  islands.  Presented  at  the 
Fourth  Conference  on  Mountain  Meteorology,  American  Meteorologi¬ 
cal  Society,  Boston,  Mass. 

FAA,  1987:  Aviation  Weather  Services,  FAA  AC-45C  (Revised),  123  pp. 

Farnham,  D. ,  and  R.  Gould,  1956:  Weather  patterns  and  local  fore¬ 
casting  at  the  Naval  ordnance  test  station,  China  Lake,  Califor¬ 

nia,  30/MS-162,  30  pp. 

Farnham,  D. ,  and  I.  Vercy,  1969:  Weather  at  the  Naval  Weapons  Center, 

China  Lake,  CA.  ,  1946-1968.  NWC  TP  4821. 

Feigenbaum,  E.  A.,  and  P.  McCorduck,  1983:  The  Fifth  Generation, 
Artificial  Intelligence  and  Japan's  Challenge  to  the  World, 
Addison-Wesley  Publishing  Co.,  New  York,  NY. 

Fingerhut  W.A.  and  P.F.  Lester,  1973:  Lower  turbulent  zones  associated 
with  mountain  lee  waves.  Paper  No.  30,  Department  of  Meteorology, 
California  State  University,  San  Jose,  California. 


Fosberg,  M.A. ,  W.E.  Marlatt,  and  L.  Krupnak,  1976:  "Estimating  Air-flow 
Patterns  over  Complex  Terrain".  USDA  Forest  Service  Rocky  Mtn. 
Forest  &  Range  Experiment  Station,  Res^ar^h  Paper 16  pp. 

Fosberg,  M.  and  M.  L.  Sestak,  1986:  KRISSY:  User’s  Guide  to  Modeling 

Three  Dimensional  Flow  in  Complex  Terrain.  PSW-92,  PSWFRES,  USFS, 
USDA,  7  pp. 

Fujita,  T. ,  1978:  Manual  of  Downburst  Identification  for  Project 

NIMROD.  SMRP  Res.  Pap._No. . 156,  University  of  Chicago,  104  pp . 

Fujita,  T.T.,  1985:  "The  Downburst".  The  University  of  Chicago,  Chicago 
Ill. 

Fujita,  T.T. ,  1986:  "DFW  Microburst".  The  University  of  Chicago, 

Chicago  Ill. 

Furman,  R.Wm. ,  and  G.L.  Wooldridge,  1984:  Observations  of  surface  flow 
over  a  three  dimensional  hill  in  a  stable  and  neutral  boundary 
layer.  Presented  at  the  Third  Conference  on  Mountain  Meteorology, 
American  Meteorological  Society,  Boston,  Mass. 

Gaffney,  J.  E.  Jr.,  and  I.  R.  Racer,  1983:  A  Learning  Interpretive 
Decision  Algorithm  for  Severe  Storm  Forecasting  Support, 
Thirteenth,  Conference  on  Seyerg  Lpcal__Stprm^^  (AMS),  17-20  Oct  83. 

GARP  Mountain  Sub-Programme,  1978:  Venice,  October  24-28,  1977,  (Copies 
of  the  report  can  be  obtained  from  Joint  Planning  Staff  for  GARP, 
c/o  World  Meteorological  Organization.  Case  Postale  No.  5,  CH- 
1211  Geneva  20,  Switzerland.) 

Gardiner,  B.A.  and  M.K.  Hill,  1984:  Use  of  an  acoustic  sounder  in  a 
mountain  environment.  Presented  at  the  Third  Conference  on 
Mountain  Meteorology,  American  Meteorological  Society,  Boston, 

Mass . 

Geiger,  R.  1965:  "The  Climate  Near  the  Ground",  Harvard  University 
Press,  Cambridge,  Mass. 

George,  J.  J.  ,  1960:  Weather  Forecasting,  for  Aeronautics.  Academic 
Press,  New  York,  673  pp. 

Giusti,  A.,  1987:  A  climatological  study  of  Alberta  Chinook  wind¬ 
storms,  associated  cross  mountain  pressure  gradients  and  upwind 
soundings.  In  Preprints  4th  Conference  on  Mountain  Meteorology, 
Seattle,  August,  1987. 

Giusti,  A.,  and  K.  P,  MacKay,  1988:  A  statistical  aid  for  the  pro¬ 
duction  of  strong  downslope  windstorms,  in  preparation. 

Greyson,  T.H. ,  G.M.  Carter,  S.  Brown,  and  A.  MacDonald,  1978:  Presented 
at  the  Conference  on  Weather  Forecasting  and  Aviation  Meteorology. 
American  Meteorological  Society,  Boston,  Mass. 


60 


Harrison,  H. ,  1966:  Mountain  wave  exposures  on  jet  routes  of  Northwest 
Airlines  and  United  Airlines.  UAL  Met.  Giro.  No.  60,  66  pp. 

Haugen,  D.  ,  (Ed),  1973:  Workshop  on  Micrometeorology ,  AMS,  Boston, 

392  pp. 

Heister  T.R.  and  W.T.  Pennell,  1980:  "The  Meteorological  Aspects  of 

Siting  Large  Wind  Turbines",  DOE  Energy  Rep.  -  Contract  EY-76-C- 
06-1830.  Battelle  Northwest  Labs,  Richland,  Wash. 

Henmi ,  T. ,  P.  Tabor,  and  T.  Lucas,  1988:  User's  Guide  to  3-d  wind  flow 
model  over  NTC  (Ft.  Irwin)  area.  5  pp  and  diagrams. 

Higgins,  H.  C.,  1963:  The  Thermal  Index,  Soaring,  21,  8-10, 

Holton,  J.R.,  1967:  The  diurnal  boundary  wind  above  sloping  terrain. 
■Tellus,  19,  199-205. 

Hovermale,  J.B. ,  1965:  "A  Non-Linear  Treatment  on  the  Problem  on 

Airflow  over  Mountains",  PhD.  Thesis,  Pennsylvania  State  Univer¬ 
sity,  University  Park. 

Hunt,  J.C.R,,  J.C.  Kaimal,  and  J.E.  Gaynor,  1985:  Some  observations  of 
turbulence  structure  in  stable  layers.  Quart.  J.  R.  Met.  Soc . , 

111,  793-815. 

Hyun,  J.M.  and  M.  Kim,  1979:  The  effect  of  nonuniform  wind  shear  on  the 
intensification  and  reflection  of  mountain  waves.  J_.  Atmos.  Sci., 
36,  2379-2384. 

Jones,  J.  R.  Mielke,  G.  Jones,  et  al.,  1970:  "Turbulence  Forecasting," 
in  Low  Altitude  Atmospheric  Turbulence  LO-LO  CAT  Phase  III. 
AFFDL-TR-70-10 ,  Volume  I,  Part  II,  60-75, 

Kessler,  E.  (Ed),  1985:  'Thunderstorm  Morphology  and  Dynamics, 
University  of  Oklahoma  Press,  Norman,  411  pp. 

Keyser,  D.,  1986:  Atmospheric  Fronts:  An  observational  perspective  in 
Mesoscale  Meteprplqgy  and.  Forecasting  (P.  Ray,  Ed.).  AMS,  Boston, 
pp  216-258. 

Klemp,  J.B.  and  D.K.  Lilly,  1975:  The  dynamics  of  wave-induced  down- 
slope  winds.  J,  Atmos.  Sci.,  32,  320-339. 

Klemp,  J.B.  and  D.K.  Lilly,  1978:  Numerical  simulation  of  hydrostatic 
mountain  waves.  J_..  Atmos  Sci .  ,  35,  78-107. 

Koch,  R.,  W.  Biggs,  P.  Hwang,  I.  Leichter,  K.  Pickerring,  E.  Sawdey, 

and  J.  Swift,  1977:  Power  Plant  Stack  Plumes  in  Complex  Terrain. 
EPAP6007-77-P20.  235  PP. 

Kuettner,  J.,  (Ed),  1986:  Scientific  Results  of  the  Alpine  Experiment 

(ALPEX) ,  Vols  I  and  II.  GARP  Pub ,  Series  No .  27.  WMO/TD  No.  108, 
710  pp  and  appendices. 


61 


Lavoie,  R.L.,  1974:  A  numerical  model  of  trade  wind  weather  over  Oahu. 
Mon.  Weath.  Rev.,  102,  630-637. 

Lee,  D.,  R.  Stull,  and  W.  Irvine,  1979:  Clear  Air  Turbulence  Fore¬ 
casting  Techniques.  USAF,  AFGWC/TN-79/001 .  73  pp. 

Lee,  J.T.,  S.  Barr,  R.E.  Lawson,  W.H.  Snyder,  and  G.L.  Marsh,  1984: 
Towing  tank  studies  of  stratified  flow  over  ridges  and  valleys. 
Presented  at  the  Third  Conference  on  Mountain  Meteorology,  Ameri¬ 
can  Meteorological  Society,  Boston,  Mass. 

Lee,  M.E.,  1988:  Low-level  turbulence  forecasting  aid.  (Unpublished) 
U.S.  Army  Atmospheric  Sciences  Laboratory,  White  Sands  Missile 
Range,  NM 

Lenschow,  P.  (ed.),  "Probing  the  Atmospheric  Boundary  Layer".  American 
Meteorological  Society,  Boston  Mass. 

Lester,  P.F.  and  W.A.  Fingerhut,  1974:  lower  turbulent  zones  associated 
with  mountain  lee  waves.  J .  Appl .  M^eor .  ,  13,  54-61. 

Lester,  P.,  R.  E,  Bach,  and  N.  Muranaka,  1988:  A  severe  Chinook 
windstorm,  in  preparation. 

Lilly,  D.K.,  1978:  A  severe  downslope  windstorm  and  aircraft 

turbulence  event  induced  by  a  mountain  wave.  JA.S,  35,  59-77. 

Lilly,  D.K,  and  J.B.  Klemp,  1979:  The  effects  of  terrain  slope  on  non¬ 
linear  hydrostatic  mountain  waves.  J.  Fluid  Mech. ,  95,  241-261. 

Lilly,  D.K.  and  J.B.  Klemp,  1980:  Comments  on  the  evolution  and 

stability  of  finite-amplitude  mountain  waves.  Part  II:  Surface 
wave  drag  and  severe  downslope  windstorms.  J.  Atmos.  Sci.,  37, 
2119-2121. 

Lilly,  D.K.  and  E.J.  Zipser,  1972:  The  Front  Range  windstorm  of  11 
January  1972.  Weatherwise,  25,  53-56. 

Lilly,  D.  K. ,  T.  Gal-Chen,  J.  Lin,  J.  Schneider,  1988:  "Phoenix  II 

Analysis  Results."  Presented  at  the  8th  Symposium  on  Turbulence 
and  Diffusion,  San  Diego,  AMS.  9  pp  and  diagrams. 

Lindsay,  C,  V.,  and  S.  Lacy,  1976:  Soaring  Meteorology  for  Fore¬ 
casters,  SSA,  Los  Angeles.  73  pp. 

Ludwig,  F.  L.,  R.  M.  Endlich,  and  K.  C.  Nitz,  1985:  Inclusion  of 
energy  considerations  in  a  mass-consistent  wind  interpolation 
scheme.  Presented  at  the  6th  Annual  EOSABL/TWI  Conference,  Las 
Cruces,  New  Mexico,  26  pp, 

Luna,  R.E.  and  H.W.  Church,  1972:  A  comparison  of  turbulence  intensity 
and  stability  ratio  measurements  to  Pasquill  stability  classes.  J. 
Appl.  Meteor.,  11,  663-669. 


Lyons,  W.A.,  1972:  The  climatology  and  prediction  of  the  Chicago  lake 
breeze.  J .  Appl . Meteor. ,  11,  1259-1270. 

Lyons,  W.A. ,  1975:  Turbulent  diffusion  and  pollutant  transport  in 
r  shoreline  environments.  {D.  Haugen,  ed.)  Lectures  on  Air  Pollution 

and  Environmental  Impact  Analysis.  American  Meteorological  Soci¬ 
ety,  Boston,  Mass.  pp.  136-208. 

*  Mass,  C.F.,  1981:  Topographically  forced  convergence  in  western 

Washington  State.  Mon.  Weath.  Rey^  109,  1335-1347. 

Mass,  C.F.,  and  M.D.  Albright,  1985:  A  Severe  windstorm  in  the  lee  of 
the  cascade  mountain  of  Washington  state.  Mon.  Weath.  Rev.,  113, 
1261-1281. 

Maher,  L.J.,  and  R.A.  Pielke,  1978:  A  test  of  an  upstream  spline 

interpolation  technique  for  the  advective  terms  in  a  numerical 
mesoscale  model.  Mon  We  a  th . ..  Re  v .  ,  106,  1758. 

Mahrt,  L.,  1985:  Vertical  structure  and  turbulence  in  the  very  stable 
boundary  layer.  J .Atmos .Sci . ,  42,  2333-2349. 

Mahrt,  L,  and  J.  Paumier,  1984:  Turbulence  over  the  coastal  range  of 
Northern  Yugoslavia.  Presented  at  the  Third  Conference  on 
Mountain  Meteorology,  American  Meteorological  Society,  Boston, 

Mass . 

Manis,  P.C.  and  B.L.  Sawford,  1979:  A  model  of  katabatic  winds.  J. 

Atmos  Sci ,  36,  619-630. 

Manis,  P.C.  and  B.L.  Sawford,  1982:  Mesoscale  observation  of  upstream 
blocking.  Quart  ..J._R^^  108 »  427-434. 

Mathews,  M.  D.,  1985:  Use  of  Computer  Technology  at  the  National 

Aviation  Weather  Advisory  Unit.  Preprint  Volume,  2d  International 
Conference  on  the  Aviation  Weather  System,  June,  1985,  Montreal, 
AMS,  Boston.  pp  285-288. 

McBean  G.A. ,  and  T.S.  Murty,  1985:  A  marked  mesoscale  pressure-induced 
wind  variation  (Research  Note).  Boundary-Layer  Met . ,  31,  209-216. 

McCarthy,  D.H.  and  J.A.  Young,  1978:  A  profile  of  the  lower  tropo¬ 
spheric  environment  of  the  New  England  coastal  front.  Presented  at 
the  Conference  on  Weather  Forecasting  and  Aviation  Meteorology . 
American  Meteorological  Society,  Boston,  Hass. 

♦  McCloskey,  J.W.,  1972:  The  effect  of  aircraft  flight  conditions  on  low 

altitude  critical  air  turbulence  (LO-LOCAT) .  Presented  at  the 
Internationa 1  Conference  on  Aerospace  and  Aeronautical  Me teorol- 
qgy,  American  Meteorological  Society,  Boston,  Mass. 

McGinley,  J. ,  1986:  "Nowcasting  Mesoscale  Phenomena"  in  Mesoscale 
Meteorology  and  Forecasting.  AMS  657-688 . 


63 


Mellor  G.L.,  and  T.  Yamada,  1974:  A  hierarchy  of  turbulence  models  for 
planetary  boundary  layers.  J.  Atmos.  Sci..  31,  1791-1806. 

Meyers,  R.E. ,  W.D.  Ohmstede,  and  R.T.  Cederwall,  1985:  A  model  of 

finite  velocity  diffusion  dynamics  from  large  eddy  simulation, 
exhibiting  both  space  and  time  memory.  Presented  at  the  Seyenth 

Symposium . p f _  Turbulence  and  Diffusion ,  Amer ican  Meteorological 

Society,  Boston,  Mass. 

Miller,  R.,  1972:  Notes  on  analysis  and  severe  storm  forecasting 
procedures,  of  AFGWC.  AWSTR-200  (Rev) ,  190  pp. 

Mitsuto,  Y. ,  and  N.  Monji,  1986:  Comparisons  of  aircraft  and  tower 

measurements  around  Tarama  island  during  the  AMTEX  '75.  J.  Climate 
Appl .  Meteor.,  25»  1946-1955. 

Monin,  A.S.  and  A.M.  Obukhov,  1953:  Dimensionless  characteristics  of 
turbulence  in  the  layer  of  atmosphere  near  the  ground.  DpklacJy 
Akademii  Nauk  .SS SR ,  93,  257-267. 

Munn,  R.  E. ,  1966:  Descriptive  Micrometeorology ,  Academic  Press,  New 
York,  245  pp. 

Nau,  D.  S.,  1983,  Expert  Computer  Systems,  IEEB_comput^er ,  16  Feb  83, 
p  63. 

Muranaka,  Neil  J.,  1988:  Application  of  a  mass  consistent  wind  model 
to  Chinook  windstorms.  M.  S.  Thesis,  Department  of  Meteorology, 
San  Jose  State  University. 

Neff,  W.D. ,  and  C.W.  King,  1985:  "Studies  of  complex  terrain  flows 

using  acoustic  sensors.".  United  States  Dept,  of  Commerce,  NOAA 
Environmental  Res.  Lab.,  Boulder,  Co.,  ASCOT  85-1. 

Nicholls,  J.  M.  ,  1973:  The  airflow  over  mountains.  Research  1958- 
1953,  WMO  TN_N0.127,  73  pp. 

NTC,  1987:  Forecaster  Handbook  for  National  Training  Center,  Bicycle 
Lake  AAF,  Ft.  Irwin,  CA.  OPR:. _ 5WS/DOJ. ,  35  pp. 

Palmen,  E. ,  and  C.  W.  Newton,  1969:  Atmospheric  Circulation  Systems, 
Academic  Press,  N.  Y.  603  pp. 

Panofsky,  H.  A.,  and  J.  Dutton,  1984:  Atmospheric  Turbulence,  Wiley- 
Interscience ,  New  York,  397  pp. 

Pasquill,  F.,  1961:  The  estimation  of  the  dispersion  of  windborne 
material.  Meteorol.  Mag,  90 (1063) ,  33-49. 

Peltier,  W.R.  and  T.L.  Clark,  1979:  The  evolution  and  stability  of 
mountain  waves.  Part  II  Surface  wave  drag  and  severe  downslope 
windstorms.  J.  At.-nos.  Sci.,  36,  1498-1529. 


64 


Peltier,  W.R.  and  T.L.  Clark,  1983:  Nonlinear  mountain  waves  in  two  and 
three  spatial  dimensions.  Quart.  J.  Roy.  Met. Spc., 109,  527-548. 

Petterssen,  S.,  1956:  Weather  Analysis  and  Forecasting,  Vol  I.  McGraw- 
Hill,  N.Y.  428  pp. 

Pettre,  P.,  1982:  On  the  problem  of  violent  valley  winds.  J.  Atmos. 
Sci.,  39,  542-554. 


Pielke,  R.A.,  1984:  Mesoscale  Meteorological  Modeling.  Academic 
Press ,  New  York,  612  pp. 

Porch,  W.M. ,  et  al.,  1988:  Tributary  contribution  to  and  interactions 

with  air  flow  in  a  deep  valley.  Presented  at  the  Fourth  Conference 
on Mountain  Meteorology ,  American  Meteorological  Society,  Boston, 
Mass . 

Queney,  P.,  G.  Corby,  N.  Gerbier,  H.  Koshmieder,  and  J.  Zierep,  1960: 
"The  Airflow  over  Mountains".  WMO  Tech.  Note  34,  135  pp. 

Racer,  1.  R. ,  and  J.  Gaffney,  Jr.,  1984:  Interpretive  Processing/ 

Expert  Systems:  An  Initiative  In  Weather  Data  Analysis  and  Fore¬ 
casting,  NWD,  Vol.  9,  No.  3,  Aug  84,  pp  31-45. 

Racer,  I.  R.,  and  J.  E.  Gaffney,  Jr,,  1986:  The  Potential  Role  of 

Artificial  Intelligence/Expert  Systems  in  the  Warning  and  Forecast 
Operations  of  the  National  Weather  Service,  ^erican_ Institute  of 
Aerospace  Sciences  24th  Aerospace  Sciences  Meeting ,  6-9  Jan  86, 

Ray,  P.S.  (ed.),  1986:  Mesoscale  Meteorology  and  Forecasting.  American 
Meteorological  Society,  Boston,  Mass.,  703  pp. 

Richardson,  L.F,,  1920:  Some  measurements  of  atmospheric  turbulence. 
Phil ,  Trans  .  Roy .  Soc  .  London ,  J^r_ A ,  221,  1-28. 

Richwien,  B.A.,  1979:  Update  on  low-level  wind  shear  forecasting,  Nat. 
^e a t h ^  Digest ,  4,  21-26. 

Sangster,  W.E. ,  1977:  "An  updated  objective  forecast  technique  for 
Colorado  downslope  winds,”  NOAA  Tech.  Memo.  NWS J3R-^  1-24. 

Scorer,  R.  ,  1978:  Environmental  Aerodynamics,  Norwood ,  England,  488  pp. 

Segal,  M.  et  al .  ,  1986:  Evaluation  of  cloud  shading  effects  on  the 

generation  and  modification  of  mesoscale  circulations.  Mon._Wea. 
Rev. ,  114,  1201-1212. 

Selvam  A.M. ,  A.S.  Ramachandra,  and  Bh.V.  Ramana,  1983:  Surface 

frictional  turbulence  as  an  agent  for  the  maintenance  and  growth 
of  large  eddies  in  the  atmospheric  planetary  boundary  layer. 
Presented  at  the  Sixth  Symposium  on  Turbulence  and  Diffusion, 
American  Meteorological  Society,  Boston,  Mass. 


65 


Sherman,  C.E. ,  1978:  A  mass-consistent  model  for  wind  fields  over 
complex  terrain.  J.  Appl.  Meteor.,  17,  312-319. 

Simpson,  R.  H. ,  1971:  The  Decision  Process  in  Hurricane  Forecasting, 
NOAA  Technical  Memorandum  NWS  SR-53.  U.S.  Dept,  of  Commerce, 
National  Oceanic  and  Atmospheric  Administration,  National  Heather 
Service,  Southern  Region,  Jan  71. 

Smith,  R.B.,  1977:  The  steepening  of  hydrostatic  mountain  waves.  J. 
Atmos.  Sci.,  34,  1634-1654. 

Smith,  R.B.,  1979:  The  influences  on  mountains  of  the  atmosphere.  Adv. 
Geophys . ,  21,  87-230. 

Smith,  R.B.,  1982:  A  differential  advection  model  of  orographic  rain. 
Mon.  Heath.  Rev. ,  110 ,  306-309. 

Smith,  R.B.,  1984:  Aerial  observations  of  the  Yugoslavian  Bora. 
Presented  at  the  Third  Conference  on  Mountain  Meteorology, 

American  Meteorological  Society,  Boston,  Mass. 

Smith,  R.B.  1985:  On  severe  downslope  winds.  J .  Atmos .  Sci.,  42, 
2597-2603. 

Smith,  S.B.  and  A.J.  Faller,  1984:  Mesoscale  flow  patterns  in  rotating 
laboratory  experiments  with  realistic  topography.  Presented  at  the 
Third  Conference  on  Mountain  Meteorology.  American  Meteorological 
Society,  Boston,  Hass. 

Stone,  G.L.  and  D.E.  Hoard,  1988:  Hind  fluctuations  on  a  nocturnal 
drainage  flow.  Presented  at  the  Fourth  Conference  on  Mountain 
Meteorology,  American  Meteorological  Society,  Boston,  Mass. 

Vinnechenko,  N. ,  N.  Pinus,  S.  Shmeter,  and  G.  Shur,  1973:  Turbulence 
in  the  Free  Atmospb«»re.  Consultants  Bureau,  New  York,  1973, 

263  pp. 

Hal ling ton,  C.  E.,  1966:  Meteorology  for  Glider  Pilots.  J.  Murray, 
London . 

Halter,  B.A.  and  J.E.  Overland,  1982:  Response  of  stratified  low  in  the 
lee  of  the  Olympic  Mountains.  Mon.  Heath.  Rev.,  110 ,  1458-1473. 

Haters,  A.,  1970:  Forecasting  Gusty  Surface  Hinds  in  the  Continental 
U.S.  AHS  TRU9.  74  pp. 

Heisman,  M. ,  and  J.  Klemp,  1986:  Characteristics  of  Isolated  Con¬ 
vective  Storms,  in  Mesoscale  Meteorology  and  Forecasting  (P.  Ray, 
Ed.)  AMS,  Boston,  pp  331-358. 

Hhiteman  C.D.,  1982:  Breakup  of  temperature  inversions  in  deep  mountain 
valleys:  Part  I  Observations.  J.  Appl.  Meteor.,  21,  270-289. 


66 


Whiteman,  C.D.,  1988:  Atmospheric  mass  budget  for  a  deep,  narrow  valley 
in  Colorado.  Presented  at  the  Fourth  Conference  on  Mountain  Meteo¬ 
rology,  American  Meteorological  Society,  Boston,  Mass. 

Whiteman,  C.D.  and  T.B.  MeKee,  1982:  Breakup  of  temperature  inversions 
in  deep  mountain  valleys:  Part  II  Thermodynamic  Model.  J.  Appl . 

‘  Meteor.,  21,  290-302. 

WMO,  1977:  Forecasting  Techniques  of  Clear  Air  Turbulence  including 
•  that  associated  with  mountain  waves.  WMO  No.  482,  31  pp. 

Wooldridge,  G.L.  and  M.M.  Orgill,  1978:  Airflow,  diffusion  and  momentum 
flux  patterns  in  a  high  mountain  valley.  Atmos .  Environ. ,  12,  803- 
808.  . . 

Wooldridge,  G.L.  and  R.Wm.  Furman,  1984:  The  use  of  Froude  numbers  to 

represent  airflow  patterns  around  an  isolated  mountain.  Presented 
at  the  Third  Conference  on  Mountain  Meteorology,  American  Meteoro¬ 
logical  Society,  Boston,  Mass. 

Wyngaard,  J.C.,  1985:  Structure  of  the  planetary  boundary  layer  and 

implications  of  its  modeling.  J.  Climate  Appl.  Meteor.,  24,  1131- 
1142. 

Wyngaard,  J.C.,  1985(b):  Large-eddy  simulation  in  small-scale 

meteorology.  Presented  at  the  Seventh  Symposium  pf_Turbjalence 
Diffusion,  American  Meteorological  Society,  Boston,  Mass. 

Yamada,  T.,  1981:  A  numerical  simulation  of  nocturnal  drainage  flow.  J. 
Meteor.  Soc.  Japan ,  59,  108-122. 

Yamada,  T.,  1983:  Simulations  of  nocturnal  drainage  flows  by  a  q21 
turbulence  closure  model.  J.  Atmos.  Sci . ,  40 ,  91-106. 

Yamada,  T. ,  and  C-Y  J.  Kao,  1986:  A  modeling  study  on  the  fair  weather 
marine  boundary  layer  of  the  GATE.  J.  Atmos.  Sci.,  43,  3186-3199. 

Yoshikado,  H. ,  1981:  Statistical  analysis  of  the  sea  breeze  pattern  in 
relation  to  general  weather  patterns.  J.  Meteor.  Soc.  Japan,  59, 
98-107. 


67 


